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Abstract 

Electron-Positron Production in Ultra-Peripheral Heavy-Ion 
Collisions with the STAR Experiment 

by 

Vladimir Borisovitch Morozov 
Doctor of Philosophy in Physics 

University of California, Berkeley 
Professor Marjorie D. Shapiro, Chair 

This thesis presents a measurement of the cross-section of the purely electromagnetic pro- 
duction of e + e~ pairs accompanied by mutual nuclear Coulomb excitation AuAu — > 
Au*Au* + e + e~, in ultra-peripheral gold-gold collisions at RHIC at the center-of-mass col- 
lision energy of y / s/v r /v = 200 GeV per nucleon. These reactions were selected by detecting 
neutron emission by the excited gold ions in the Zero Degree Calorimeters. The charged 
tracks in the e + e~ events were reconstructed with the STAR Time Projection Chamber. 

The detector acceptance limits the kinematical range of the observed e + e~ pairs; 
therefore the measured cross-section is extrapolated to 4ir with the use of Monte Carlo 
simulations. We have developed a Monte Carlo simulation for ultra-peripheral e + e~ pro- 
duction at RHIC based on the Equivalent Photon Approximation, the lowest-order QED 
e + e~ production cross-section by two real photons and the assumption that the mutual nu- 
clear excitations and the e + e~ production are independent (EPA model), various kinematic 
regions. 

We compare our experimental results to two models: the EPA model and a model 
based on full QED calculation of the e + e~ production, taking the photon virtuality into 
account. The measured differential cross-section da/dMi nv (Mi nv - e + e~ invariant mass) 
agrees well with both theoretical models. The measured differential cross-section da/dp^ 1 
(jptot _ e + e - total transverse momentum) favors the full QED calculation over the EPA 
model. 
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Professor Marjorie D. Shapiro 
Dissertation Committee Chair 



"The central and peripheral collisions of relativistic heavy ions may be com- 
pared to the case of two potential lovers walking on the same side of the street, but 
in opposite directions. If they do not care, they can collide frontally... It could be 
a good opportunity for the beginning of strong interactions between them. ... On 
the other hand, if they pass far from each other, they can still exchange glances (just 
electromagnetic interaction!), which can even lead to their excitation. ...the effects 
of these peripheral collisions are sometimes more interesting than the violent frontal 
ones." 



G. Baur and C. A. Bertulani 
Physics Reports 163, 299, (1989) 
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Chapter 1 

Introduction 

This thesis presents a study of the production of e + e~ pairs in ultra-peripheral 
AuAu collisions at the Relativistic Heavy Ion Collider (RHIC) observed with the STAR 
detector (Solenoidal Tracker at RHIC). While most of the RHIC physics program is con- 
cerned with the hadronic AuAu interactions, we focus on the ultra-peripheral interactions, 
where the Au ions interact only via the long-range forces. The e + e~ pairs are a product of 
purely electromagnetic interaction of the virtual photon fields emitted by the Au ions. The 
process of two photons converting into an e + e~ pair has been studied experimentally with 
a good precision in quantum electrodynamics (QED) for pair energies up to 100 GeV. 

The study of this process at RHIC, however, opens a door for several new interest- 
ing features. Since the fields from the proton constituents of Au ions (with charge Z = 79) 
add coherently, the electromagnetic field strength scales as Z 2 and the interaction rate as 
Z 4 . Therefore, RHIC will produce a very large two-photon interaction rate. Additionally, 
the field coupling constant is scaled up by the charge of Au ions, Za ~ 0.6, and the photon- 
photon interactions enter a strong interaction regime. We wish to measure the cross-section 
of the e + e~ production in this regime. 

The STAR detector is optimized for detecting charged particles at mid-rapidity 
and in the transverse momentum range from about 100 MeV/c up to several GeV/c. The 
fraction of e + e~ pairs with both electron and positron tracks reaching the STAR Central 
Trigger Barrel (located at mid-rapidity) is close to zero; therefore triggering on the exclusive 
reaction AuAu — ► AuAu + e + e~ in STAR is very difficult. We chose to focus our attention 
on a closely related reaction - electromagnetic production of an e + e~ pair with mutual 
nuclear Coulomb excitation of the colliding ions: AuAu — ► Au* Au* + e + e~ . The neutrons 
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Au Au 




Figure 1.1: Electron-positron production with a mutual Coulomb excitation. Thee + e~ pair 
production (right of the dashed line) is independent of the simultaneous nuclear excitation 
(left of the dashed line). 

emitted when the nuclear excitations decay provide a tag on this reaction in the forward 
neutron calorimeters. 

Figure 11.11 gives a schematic view of this process. Each Au ion emits a field of 
virtual photons, but this process doesn't disrupt the emitting ion. The virtual photons can 
produce an e + e~ pair in collisions with the photons from the ion in the opposite beam. 
Additionally, photons can cause an excitation of the Au ions in the opposite beam. We 
trigger on the AuAu collisions in which both Au ions are excited via a photon exchange. 

Recent theoretical and experimental studies of mutual Coulomb excitation in 
heavy-ion collisions suggest that mutual excitation is independent of other reactions in 
the collision. However, Coulomb excitation has not been observed in coincidence with 
electromagnetic e + e~ pair production in any previous heavy- ion collision experiments. 

The study of the reaction AuAu — ► Au* Au* + e + e~ in STAR thus gives us a 
chance to examine two important physics questions. First is the test of QED at high fields 
- are there indications that at the available energies more than two photons are involved in 
an e + e~ pair creation? The second question is whether pair production is independent of 
photonuclear excitation. 
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Chapter 2 

Physics of Electron-Positron Pair 
Production 

Electron-positron pair production is a purely electromagnetic process; therefore 
we start with a brief overview of quantum electrodynamics (QED). We then discuss the 
special features of e + e~ pair production in relativistic heavy-ion collisions. We present 
the lowest-order QED framework for computing the e + e" production cross-section, and 
discuss its higher-order extensions. Finally we present a model for the computation of the 
AuAu — > Au*Au* + e + e~ cross-section. 

2.1 Overview of Quantum Electrodynamics 

Quantum electrodynamics is one of the simplest Abelian gauge quantum field 
theories of nature. The theory incorporates the ideas of Maxwell's electromagnetism with 
quantum mechanics, such as quantization and the spin-1/2 nature of electrons and positrons. 
The key results of the theory can be summarized as follows: 

- The electromagnetic field consists of photons - the massless spin-1 corpuscles of energy 
of the field. 

- The interaction of electrons and positrons with the field (and each other) can be 
represented as the interaction between the electrons/positrons and the photons. 

- The electrons/positrons and photons can be born or annihilated in interactions with 
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each other. 
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In accordance with this idea, the QED Lagrangian consists of three terms: a term 
from relativistic quantum mechanics describing free electrons/positrons, a term describing 
the free Maxwell field, and a term responsible for the interaction between electrons/positrons 
and the field: 

^QED = ^electron + -^Maxwell + -^Interact = 1p (id - 77t e ) ijj - - {F^f - eip^lpAfj, (2.1) 

where ip is the Dirac wave- function of electrons/positrons, is the electromagnetic vector 
potential, F^ v = d^A u — d u A^ is the electromagnetic field tensor, and e and m e are the 
charge and the mass of the electron. We use the natural system of units where c = 1 and 
h = l. 

The Lagrangian (|2.1|) considers only one fermion family - electrons and positrons. 
Amazingly, this very simple Lagrangian can account for nearly all observed phenomena from 
macroscopic scales down to 10~ 13 cm. At distance scales below 10~ 13 cm (or, equivalently, 
for interaction energies above 100 GeV) phenomena explained by other theories, such as 
electro- weak theory or quantum chromo dynamics become significant. 

One of the central applications of the QED Lagrangian is to find the transition 
amplitude for some initial state which, generally speaking, consists of particles 1, 2, etc. 
with momenta p\, p 2 , etc. at time to (denoted | Initial)), to make a transition into a final 
state consisting of particles A, B, etc. with momenta pa, Pb, etc. at time t\ (denoted 
(Final|). According to the rules of quantum mechanics the transition amplitude is: 

(Final | Initial) = (px,p 2 ...\ exp dt J d 3 x ■ L Interact (t, x)| \p A ,PB-) (2.2) 

where the expression sandwiched between the initial and final momenta is the usual time- 
evolution operator, based on the interaction part of the QED Lagrangian ()2.1|) . The com- 
putation of cross-sections for all QED processes is based on variations of formula H2.2JI . 1 

2.1.1 Perturbative Expansions in Quantum Electrodynamics 

The analytical computation of the transition amplitude (|2.2j) is impossible for most 
non-trivial initial/final state combinations. Fortunately, in 1949 Fevnman|44| proposed that 
the exponential in the time-evolution operator can be expanded as: 



exp j-i J dtd 3 x ■ Li ntcj:act (t,x)^ 



1 The formula given here is for conceptual demonstration only, for details see the standard texts 




Figure 2.1: Lowest order Feynman diagrams for: a) Bhabha scattering; b) Compton scat- 
tering; c) pair annihilation; d) pair creation. Topologically identical diagrams not presented 
in b), c) and d). 

J dtd 3 x-Li Qteract (t,x)+(-i) 2 jj dt 1 d 3 x 1 dt 2 d 3 X2-Li nteract (t 1 ,x 1 )-Li QteTSuCt (t2,X2)+..= 



= l+(-z)v / 47ra J dtd^xip^ipA^+Anai-i) 2 J J dt 1 d 3 xidt 2 d 3 x 2 ijJ 1 ^ipiA lfi • ip^^A^... 

(2.3) 

where a = e 2 /47r ~ 1/137 is the fine structure constant. The fine structure constant, often 
referred to as the coupling constant, represents the coupling strength of the electromagnetic 
field to the electron charge. Effectively Equation (|2.3fl is a Taylor expansion in the orders 
of a 1 ' 2 . 

For the non-identical initial and final states, the lowest-order non-zero matrix ele- 
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Figure 2.2: Pair kinematics: a) longitudinal projection in the lab frame, b) transverse 
projection in the lab frame, c) longitudinal projection in the center of mass frame of the 77 
pair. 

ment in this perturbative expansion is of the order a|58j. Figure I2~T1 shows the lowest-order 
QED Feynman diagrams for: e + e~ — ► e + e~ (Bhabha scattering), e~7 — ► e~7 (Comp- 
ton scattering), e + e~ — > 77 (pair annihilation) and the reverse process 77 — > e + e~ (pair 
creation). 2 

The lowest order is a good approximation to the total amplitude if the energies 
of the particles involved in the reaction are small, so that each n-th term of the power 
expansion is dominated by the small coefficient a™/ 2 . However, the higher order corrections 
can be measured in high-precision experiments, and provide further confirmation to the 
QED 54 . In fact, to this day QED is the most stringently tested and the most successful of 
all physical theories, with agreement between the theory and the data up to eight significant 
digits. 

2.1.2 Kinematical Properties of 77 — >e + e~ Process 

Let's consider the properties of two-photon annihilation into an e + e~ pair closely. 
We assume that the two colliding photons have very large momenta pi~ and p2-y of opposite 
signs along the axis z and very small momenta in the plane transverse to the z axis. 3 Figure 
12.21 shows the kinematical variables used in the description of this reaction. 

The six variables (Px y x and p 2 i2/)Z ) describing the individual track kinematics 

uniquely define the six variables describing the pair kinematics: pair invariant mass (Mi nv ), 

2 The processes involving only photons in the initial state (e.g. Figure ETT1 (d)) are referred to as 'two- 
photon physics'. 

3 This will always be the case for the photons produced in the heavy-ion collisions, as we explain later. 
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pair rapidity (Y), absolute value of the pair total transverse momentum (p*£ ), the azimuthal 
angle of the total transverse momentum (<p), 4 the angle between the total momentum and 
the momentum of the negative track in transverse plane (ip), 4 and the polar angle of the 
negative track in the center-of-mass frame (6'). The formulae for computing these variables 
from the two tracks' momenta are provided in Appendix lAl 

In the center-of-mass frame, the kinematics of the reaction is defined by the angle 
between the track momenta and the direction of the photons 9* (Figure T2. 21 c). This angle is 
the same for the electron and positron tracks. Determination of the 9* generally requires the 
knowledge of photon momenta (p^x^P^y^P^z)] however, if the photon transverse momenta 
are approximately zero, 9* can be approximated by 9'. We use a Monte Carlo simulation 
in Chapter |3] to examine the accuracy of this approximation. In the lab frame the angles 
between the track momenta and the z axis are different {9 e + and 9 e -). Traditionally, 9 e + 
(or 9 e - ) is expressed via pseudorapidity, defined as 

v track = _ Q g ]n (l etrack \ \ (2.4) 

The computation of the cross-section for the e + e~ pair production (in the pair's 
center-of-mass frame) in the lowest QED order is performed according to the Feynman 
rules for the diagram in Figure l2~Tl (d)^Hl- F° r t w ° colliding circularly polarized non- 
virtual photons the angular distribution of the cross-section in the center-of-mass frame 
is: 



da + e - 2 2 Ms - 4m 2 ) 

— Ion a • • 



dcos9* V s 



,3 



yfs — \fs — 4m 2 cos 9* yfs + \fs~— 4m 2 cos I 



+ \fs — 4m 2 cos 9* — Vs - 4m 2 cos 9* 

(2.5) 

where s = Mf nv is the squared invariant mass of the system, m is the mass of the electron, 
and a is the fine structure constant. Appendix^]discusses the effects of photon polarization, 
which are negligible in the rangle of energies available at RHIC. 

For the large 77 system invariant masses (m 2 <C s) the right side of Equation (|2,5|) 
can be approximated as da/dcos9* ~ (1 + cos 2 9*)/ sin 2 9*. Thus the electron/positron 
production is peaked in the direction of the incoming photons. 

Integrating the angular cross-section (|2,5|) one obtains the full cross-section to 
4 For non-degenerate p 1 " 1 / cases. 



9 



produce an e + e pair: 

5,1/ 4m 2 8m 4 \, ( J~s r~s \ , (\ 4m 2 \ / 4m 2 
^ e+e -( S ) = 8- 2 -(l + ^-^jln^ + ^-lj + 4vra 2 (- + —) — 

(2.6) 

As a general property of all 77 processes, the cross-section depends only on y/s. In their 
center-of-mass frame the two photons have equal energies and therefore only one variable 
defines the reaction. 

The cross-section of the e + e~ pair production has a peak at s ~ (2.5m) 2 , and 
drops asymptotically as ln(s)/s for y/s 3> m. The differential cross-section for e + e~ pair 
production is maximal at the invariant masses close to the mass of the electron and in the 
forward region. 

2.2 QED in Relativistic Heavy-Ion Collisions 

To understand e + e~ with nuclear excitation, it is useful to first consider exclu- 
sive e + e~ production. We will first describe the fundamental approach to the two-photon 
processes in the ultra-peripheral relativistic heavy-ion collisions - the Equivalent Photon 
Approximation. This approximation uses the lowest-order term of perturbative QED. Next 
we will discuss the corrections to the cross-section from the higher-order terms and from 
non-perturbative calculations. Finally, we will discuss the previous experimental measure- 
ments of e + e~ pair production at e + e~ colliders and heavy- ion colliders. 

2.2.1 Pair Production in Ultra-Peripheral Heavy-Ion Collisions in the 
Lowest-order QED Approximation 

The basic framework for studying the two-photon physics in the high-energy col- 
lisions was first introduced by Fermi and then developed in detail by Landau and 
LifshitzjSH]. Their idea was that high-energy projectiles emitting the electromagnetic fields 
retain almost all of their initial momenta in the ultra-peripheral interaction, and thus can be 
assumed to follow classical straight-line trajectories with constant velocities. Therefore, the 
ultra-peripheral process can be thought of as a convolution of three processes: an emission 
of a photon by one of the colliding particles, photon emission by the other particle, and 
annihilation of two photons into a final state (labeled X)[17j. The two photons' densities 
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Figure 2.3: Equivalent Photon Approximation. Coupling of the fields to the Au ions is 
assumed to be classical, b is the impact parameter - the distance between the colliding 
ions. The exchange photons have four-momenta and q%. 



can be combined into a new quantity: the two-photon luminosity L 77 per ion pair. Using 
the two-photon luminosity, the cross-section of the two-photon production of the state X 
in the ultra-peripheral collision of two heavy ions (AA) becomes: 



d(JAA-*AA+X _ dL 17 



a 7 ^ x (M inv ,Y) (2.7) 



dM inv dY dM inv dY 
where M^ nv and Y are the invariant mass and rapidity of the system. 

Figure l2~Tfl demonstrates the idea graphically for the case of e + e~ pair production. 
The exchange photons have four-momenta q^ = (uj,p±,p z = ui//3), where lo is the energy of 
the exchange photon and (3 is the velocity of the nucleus. The virtuality of such photons is 
Q 2 = —q 2 = uj 2 /j 2 + p\, and 7 = — (3 2 ) is the Lorentz boost of the nucleus in the 

lab frame. 



Equivalent Photon Approximation 

Ignoring the photon virtuality, the photon flux accompanying a heavy nucleus is 

given by the Weizsacker Williams approximationjSJ. This approximation considers the ions 

as point-like sources of electric field, 5 and the field distribution in the lab frame is taken as 

a field of the charge Ze moving with the speed . For | j3 | ~ 1 the E and B fields are 

mostly in the plane transverse to the ion's velocity, with vector E pointing radially away 

from the ion and vector B perpendicular to E . Thus at each point in the transverse plane, 

5 Since we consider only the fields outside of the ions, the exact details of the charge distribution inside 
the ions are unimportant for the calculation of the photon energy spectrum. 
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the (E , B) field combination can be approximated by a plane electromagnetic wave, prop- 
agating in the direction of (5 . The equivalent photon spectrum (energy per unit frequency 
interval per unit area) as a function of the distance to the ion b is|51j: 

d 2 b V y ^ \7/ P lP 
where K\ is the modified Bessel function of order one. The Bessel function becomes very 
small for x > 1. Since we are considering fields outside of the nuclei (b > R), this means that 
the photon energies are limited to u < j/R, where R is the nuclear radius. The physical 
interpretation of this cutoff is that in order for the electromagnetic field to couple coherently 
to the ions, the photon wavelength in the rest frame of the ions should be greater than the 
size of the ion: A > R. Thus, in the nucleus rest frame coherence limits the transverse and 
longitudinal momenta of exchange photons to p± < 1/R ( ~ 30 MeV/c for Au), and in 
the lab frame (with the Lorentz boost in the longitudinal direction) the photon energy is 
limited to u < j/R sa 3 GeV > p ± . 

To calculate an overall two-photon density F 77 we need to integrate expression 
()2.8|) over all space and over all impact parameters. This introduces one important com- 
plication. When the nuclei physically collide (b < 2R), then the hadronic interactions 
will completely overshadow the electromagnetic ones. Therefore, in calculating the usable 
two-photon density, this overlap region must be excluded |13j: 



oo oo 

F 77 = J d 2 h J d 2 ^A^(?!,u;i) -^-(Ta.wa) ■&(\h-h\ -2R) (2.9) 

RAu RAu 

the ^-function above reflects an exclusion of the nuclear overlap region from the density 
calculation. 

We can convert the two photon energies wi , 102 into the invariant mass and the 
rapidity of the 77 pair in the lab frame: Mi nv = 2^/cjiu>2 and Y = l/21og(cJi/u;2), an d the 
two-photon differential luminosity can be expressed as a function of these two variables: 6 



(M m „,y) = ^F 77 f^exp(y),^exp(-Y)) (2.10) 



dM im ,dY v ' M ir ,,> "V 2 " v " 2 



3 For the derivation, see Appendix 1X1 
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Figure 2.4: Two photon differential luminosities vs. 77 pair invariant mass M 77 at RHIC 
for gold (solid), iodine (dashed) and oxygen (dot-dashed) beams p2*] . 

Klein and Nystrand in |H2] used the above method to calculate total two-photon 

differential luminosity dLtfy / 'dMi nv = L^adronic ' / ^M- ^dY ^ a * RHIC for several ion 

—00 

species for the photon center-of-mass energies up to 6 GeV. For gold beams 'jAu = 108 
and LauAu = 2.0 • 10 26 cm _2 s _1 design RHIC Lorentz boost and luminosity are assumed 
( 7/ = 104 and L n = 2.7 • 10 27 cni'V 1 for iodine, 7o = 135 and L o = 9.8 • 10 28 cm^s" 1 
for oxygen). Figure 12.41 shows the results of the calculation. The highest two-photon 
luminosity is obtained in I+I interactions. The lower Z is compensated by higher nuclear 
luminosity and smaller nuclear radius. The two-photon luminosity drops off quickly with 
increasing two-photon center-of-mass energy, which contributes to the rapidly decreasing 
total cross-section (7AuAu^AuAu+x as a function of the state X invariant mass. 

Applicability of the Lowest-order Approximation for e + e Production at RHIC 

The use of the lowest-order QED approximation and the EPA explicitly assumes 
that the e + e~ pair is produced by exactly two real photons. Photon virtuality at RHIC 
is typically on the order of Q 2 ~ l/R Au = (30) 2 MeV 2 . This is much greater than the 
electron mass squared (m 2 = (0.51) 2 MeV 2 ), and the photon virtuality cannot be neglected 
for some regions of electron-positron phase space. Additionally, the field coupling to the 
Au ions is very strong - Za ~ 0.6 at RHIC - and taking only the lowest-order term in 
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the perturbative expansion Q2.3|) might not be accurate. Lastly, the Compton wavelength 
of the electron = l/m e = 386 fm is much larger than the typical impact parameter at 
RHIC for exclusive pair production (~ 100 fm), and the e + e~ pair production is poorly 
localized. The following subsection gives an overview of the theoretical papers dealing with 
these complications, and presents a prediction of the e + e~ production cross-section in the 
lowest QED approximation and corrections from higher-order terms. 

2.2.2 Electron-positron Production Cross— section at RHIC 

The total e + e~ pair production cross-section in ultra-peripheral heavy- ion colli- 
sions can be computed as a convolution of the two-photon luminosity and the cross-section 
for 77 — ► e + e~ . To avoid problems with photon virtuality, it is customary to set the lower 
integration limits in Equation (|2.9j) equal to instead of Rau- It has been shown that the 
total e + e~ pair production cross-section is dominated by large impact parameters, there- 
fore cutting out the region b < \ e c does not affect the cross-section significantly. Alscher et 
al. in [2j compute the resulting total cross-section for the exclusive e + e~ pair production 
at RHIC at yfsNN = 200 GeV to be 33 kb. However, most of the particles are produced at 
low invariant masses (below 10 MeV) and into the very forward direction, so the fraction 
of the cross-section visible to the detectors is very small. We will estimate the observable 
cross-section in the next chapter. 

At impact parameters below this approach breaks down for the e + e~ pairs. 
Formally, we can write the probability to produce an e + e~ pair if colliding ions are at 
impact parameter b as: 

oo oo 

f f dnjut, b) dn(u 2 , b) du x dw 2 , . , . 

P{b)= J J ~Tb ^rr(^i^) (2-H) 

^min ^min 

where the photon densities n(u,b) are from ()2.8|) . and the cross-section a 77 is from (|2.6|) . 
Since the two-photon cross-section is highly peaked at ~ 2m e , we can approximate a™(ui, UJ2) ~ 
<7 77 (m e ,m e ) ~ a 2 /ml and assume u> min = m e . With this modification, Q2.11|) becomes: 

00 2 00 2 

m K / ^, (^) k\ (^) ^ / ^ (^) ki (^)^ (, 12 , 

00 

Using the identity J xK\ (xa) dx = -^\n (^) , where 5 « 0.681 is a constant, related 



14 



to Euler constant, P(b) is proportional to: 



P (6)a(Zo) <_L^(i2_) 



(2.13) 



A more careful calculation presented, for example, in yields: 




2m e b 



■yS 



) 



(2.14) 



where 7 = 27L ab — 1 is the Lorentz factor in the target frame. Equation ((2.14(1 shows that 
at RHIC energies (7Lab = 108) and for impact parameters less than the \ e c , this probability 
exceeds 1 (unitarity violation) , which clearly demonstrates a breakdown of the lowest order 
QED approach for e + e~ production at RHIC. 

Further analysis (11 j showed that production of multiple e + e~ pairs (e.g. Figure 
12.51 d)) is significant in the energy range where EPA prediction violates unitarity. Including 
multiple pairs in the cross-section computation restores unitarity. The pair production 
probability is then distributed according to a Poisson distribution: 



where lowest order contribution P^ 1 ' (b) is the mean of the Poisson distribution. Multiple 
pair production is dominant over the single pair production for b < A^. at RHIC energies. 
For AuAu collisions at ^snn = 200 GeV at RHIC, each ultra-peripheral e + e~ pair in the 
invariant mass range Mj nt) ~ 10 4- 20 MeV is expected to be accompanied on the average 
by one more e + e~ pair. However, the second pair should be typically low-energy, and not 
visible to the detectors [3]. 



a number of complications arise. There are many higher-order Feynman diagrams for this 
reaction. Two of such possible diagrams including an exchange of three and four photons are 
presented in Figure l2~5l (b^ and (c). Adding an extra photon to the higher-order corrections 
diagram suppresses the contribution by Za ~ 0.6 which is not significantly smaller than 
the previous lowest-order contribution. 



couples to the Au ion by more than one photon, make an especially large contribution. This 
correction is historically called 'Coulomb correction', since it represents a re-scattering of 
the electron or positron in the field of the nucleus. If observable, the Coulomb corrections 



P N {b) 




(2.15) 



If one wishes to compute higher-order QED contributions to the e + e cross-section 



Diagrams like the one in Figure l2~5l (b), where either the electron or the positron 
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a) Au + b) Au + c) Au+ d) Au + 




Au + Au + Au + Au 



Figure 2.5: a) Lowest order Feynman diagram for the e + e~ production b) next to lowest 
order three-photon exchange diagram c) four-photon exchange diagram, d) double e + e~ 
pair production dominant contribution. 



should make the distributions of da e + e - / dp e + and da e + e - / 'dp e - (p is the momentum of the 
electron or positron) to be different, since the electron and positron will scatter differently 
in the field of the positively charged Au ions. 

Authors of ^2] approached the computation of the Coulomb corrections in the 
collisions of charges Z\e and Z^e using the results of the Bethe-Heitler process 7 + Z —* 
e + e~ + Z where higher-order effects are well-known. Such a treatment requires an ad hoc 
symmetrization with respect to Z\ and Z2 and the results of this analysis are being disputed 
by several authors [BTI] . 

Some authors abandoned the Feynman perturbative expansion in Equation Q2.3JI 
and solved the Dirac's equation for the e + e~ in the electromagnetic field of Au ions non- 
perturbatively. In this approach the authors work in the light-cone coordinates where the 
effect of electromagnetic fields has a form of a delta-function. Thus the Green's function for 
the exact wave function of e + e~ at the interaction point is found. The transition amplitude 
is constructed from the Green's functions. The results were found to match the perturbative 
calculations with Coulomb correction [221 1551 128j . 

The agreement of the non-perturbative calculations for a single pair production 
with the perturbative calculations is somewhat surprising, given the large coupling constant 
Za. However, it was also observed that for multiple pair production, the non-perturbative 
result was smaller than the perturbative result f55]. Since multiple pair production is natu- 
rally a higher order process, it's not surprising that a difference appears. 
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7=E/M 

Figure 2.6: Comparison of e + e~ production cross-section in the lowest QED order ('Born 
cross-section') and cross-section with Coulomb corrections as a function of beam Lorentz 
factor 

Comparing Lowest-order and Non-perturbative Calculations for RHIC 

Using the perturbation theory to all orders the authors of [Slj find that Coulomb 
corrections to the lowest-order total e + e~ cross-section are negative and equal to ~ —25% 
for Au ions at ^smn = 200 GeV. Figure l2~6l compares the lowest-order QED cross-section 
and the cross-section with Coulomb correction as a function of beam energy. 

The total cross-section of the e + e~ pair production is dominated by the low 
invariant mass pairs (Mi nv ~ 2m e ). However, STAR detector can only observe e + e~ pairs 
above a certain minimal invariant mass cutoff (M™°_ S> 2m e ), and the magnitude of 
Coulomb corrections in the STAR observable range is unclear. Qualitatively, we expect 
that the Coulomb corrections become less significant for such pairs. Compared to the 
lowest-order Feynman diagram, the Coulomb correction diagrams include at least one extra 
photon propagator. The photon propagator suppreses the diagram by the factor ~ l/g 2 j5H]. 
Since high invariant mass pairs are created by high energy exchange photons, the virtuality 
of the exchange photons is typically higher for the high-Mj n?; pairs than for the low-Mj„„ 
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pairs (virtuality Q 2 = oj 2 /^ 2 + £>j_)- Therefore, adding an extra photon to the Feynman 
diagram suppresses higher-order contributions more strongly for the high-Mj ni , pairs than 
for the low- Mi nv pairs. However, the exact computation of Coulomb corrections for high- 
Mi nv pairs is not yet available 6 . 

2.2.3 Transverse Momenta in Equivalent Photon Approximation 

The equivalent photon approximation assumes zero transverse momenta for the 
equivalent photons. To estimate the photon transverse momentum spectra authors of ^3] 
consider a coupling of a photon with 4-momentum to the ion charge with an extended 
form factor F(t). The number of photons with energy oj and transverse momentum p± is 
given by: 

ocp ± — — (2.16) 



cLudp± \ t 

where t = —q 2 = u 2 j^ 2 + p\ and F(t) is a Fourier transform of the nuclear charge density 
p(r). The nuclear charge density is assumed to have a Woods-Saxon distribution: 

P(r) oc ) s ( 2 - 17 ) 

1 + exp ' r -=^ 



where nuclear radius is calculated from Rnuc = VqA 1 ^ 3 with r Q = 1.16(1. — 1.16j4~ 2 / 3 ) fm, 
and a constant skin-thickness of a = 0.535 fm is used. These parametrizations are obtained 
from electron-nucleus scattering datajlUj. 

For a heavy nucleus p(r) can be well approximated by a convolution of a hard- 
sphere and a Yukawa potential. The Fourier transform of F(t) is then a product of separate 
Fourier-transforms 19 : 

F(t) oc t^l 2 (sin (t^RNuc) - t 1/2 i?jv uc cos (t^ 2 R Nve )) ■ ( - —A where c = 0.7 



1 + tc 2 

(2.18) 

The distribution of the photons as a function of the photon energy lo and transverse 
momentum p± is then: 

dN ( u , p± ) sin ( y/tR Nuc ) - y/tR Nuc cos ( VtR Nuc ) p\ 1 



{(ojM 2 +P 2 ± ) 3/2 {(lo/ 7 ) 2 +pI) A [i + ^ h )2 +p 2_f c 2^ 

(2.19) 



where 7 = 108 is the Lorentz factor in the lab frame. Equation (|2.19|) shows that for a 
photon of energy u the perpendicular momentum cannot be much higher than lo/j, or 
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Figure 2.7: Left: L3 observed invariant mass spectrum of e + e pairs. Right: L3 measured 
differential cross-section do jd cos 6* [31] . 



w/108. The transverse momenta of the photons are negligible compared to the photon 
energies, and the use of p± = approximation in the calculation of the e + e~ production 
cross-section is justified. 

In general, the low transverse momentum of the equivalent photons and the low 
total transverse momentum of the produced final state is one of the defining characteristics 
of ultra-peripheral interactions at heavy-ion accelerators. We will use this property for 
separating the true signal from the background in our analysis. 



2.2.4 Previous Measurements 

The subject of two-photon physics first received wide attention at the Kiev Confer- 
ence in 1970, with the reports of Brodsky JH] and BalakinpJ]- In 1971 the first observation of 
'electroproduction' reaction e + e~ — > e + e~ + e + e~ was made at the colliding beam machine 
VEPP-2 in Novosibirsk |42j . This experiment measured about 100 of e + e~ pairs produced 
at the colliding beams energies of ~ 510 MeV, and the cross-section in the acceptance region 
was found to be 20.0 mb. 

Since then a large number of experiments measured two-photon e + e~ produc- 
tion at electron-positron colliders. L3 experiment at LEP collider observed the untagged 
reactions e + e~ — > e + e~ + e + e~ at = 91 GeV. The total available integrated luminosity 
was ~ 52 pb -1 [Hl], which yielded good event statistics. Figure l2~Tl compares Monte Carlo 
predictions to the observed e + e~ pair invariant mass (W 77 ) spectrum and the differential 
cross-section da jd cos 6*. Due to detector acceptance the observable kinematical range was 
limited to (44° < 8* < 136°) and W-y-y > 500 MeV. The total observable cross-section was 
measured to be (Jdata = 2.56 ± 0.01(stat.) ± 0.05(syst.) nb, which is in excellent agreement 
7 9* is defined as lepton polar angle in the center-of-mass frame (equivalent to angle 6' in Figure |2|2j 
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Figure 2.8: CERES measured differential cross-sections for e + e pairs. Left: da /dm, 
center: da/d(j} op , right: da/dp±. Dots - data, solid line - Monte Carlo |3*5] . 



with the lowest-order QED Monte Carlo. 

The first measurement of the electromagnetic e + e~ production in heavy-ion col- 
lisions was made by the CERES /NA45 Collaboration 38 . CERES/NA45 is a fixed-target 
experiment dedicated to the measurement of e + e~ production in proton-proton, proton- 
nucleus and nucleus-nucleus collisions. For the ultra-peripheral e + e~ production studies 
sulphur beams incident on lead target were used ( 16 S32 and 82 Pb207)> with the beam energy 
of 200 GeV per nucleon. Only events with no hadronic interactions taking place were ac- 
cepted for the study. Figure 12.81 shows the differential cross-section as a function of the pair 
invariant mass m ee , pair transverse momentum p± and the pair azimuthal opening angle 
(j) op . The detector acceptance limited the kinematic variables to 10 MeV < m ee < 100 MeV, 
E e ± > 25 MeV, 141 mrad < 6> e ± < 260 mrad and 90° < 4> op < 180°. 8 The data agreement 
with the lowest-order QED Monte Carlo is good for the invariant mass distribution. The 
p± distribution shows a slight enhancement for p± > 20 MeV/c, which correlates with a 
slight excess at azimuthal opening angles 4> op < 120°. This disagreement is attributed to 
statistical uncertainties and possibly imperfect background subtraction. The total cross- 
section was found to be a = 13.9 ± 3.1(stat.) + 3.5/— 1.9(syst.) mb, which is in very good 
agreement with the lowest-order QED prediction. 

^Etrack is track energy, and Qtrach is the track polar angle in the lab frame. 
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Figure 2.9: e + e~ production with mutual nuclear excitation. We assume that the pair 
production (right of the dashed line) is independent of the nuclear excitation (left of the 
dashed line). 

2.3 Electron-Positron Pairs with Mutual Coulomb Excita- 
tion in Relativistic Heavy-Ion Collisions 

The gold ions moving in the accelerator beams are surrounded by the flux of virtual 
photons. In addition to producing an e + e~ pair these photons are capable of exciting gold 
ions from the opposite beam. The diagram in Figure I2~U1 shows a production of an e + e~ pair 
accompanied by an exchange of two photons between the Au ions, leading to the mutual 
Coulomb excitation of the ions. The leading mode of the excitation is excitation into the 
state of the Giant Dipole Resonance (GDR) j^j. This collective excitation usually decays by 
a single neutron emission. Upon emission the neutrons move in the longitudinal direction 
with approximately the same momentum as the beam. 



The diagram in Figure l2~§l corresponding to the e + e pair production by two 



photons plus an exchange of two photons, is the dominant lowest-order diagram for e + e~ 
pair production with Coulomb breakup. In 1955 S. N. Gupta demonstrated that emission 
of photons involved in e + e~ pair production is independent of the emission of both pho- 
tons which are involved in the excitation of the ions |47|. This means that the probability 
of the e + e~ production at a given impact parameter b in the lowest order is indepen- 
dent of the probability of the excitation and breakup of the ions. This property is called 
factorization |49l @J. If factorization holds, the total cross-section for the ultra-peripheral 
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e + e production with simultaneous nuclear breakup is: 



a e+e~ with mutual breakup — J d 2 b P e + e -( b )PxnXn( b )P no hadronic( b ) (2.20) 

where P e + e - is the probability for an e + e~ pair production, P no hadronic is a probability 
of no hadronic interaction happening between the nuclei, and Px-nXn is a probability of a 
simultaneous nuclear excitation with breakup. 9 

The quantity P no hadronic (b) modifies the ^-function cut-off previously used in Equa- 
tion (|2.9f) . While the theta- function represents a hard-sphere nuclear model, the compu- 
tation of P no hadronic assumes a nuclear form-factor with smooth edges. Using a Glauber 
Model we get: 

-Pno hadronic(fe) = exp ^-a nn j Ta{j)T b { b - ~r t )d 2 ~r'^j (2.21) 

where a nn = 52 mb is the hadronic nucleon-nucleon cross-section (at ^snn = 200 GeV) and 
T(r) is a nuclear thickness function. The nuclear thickness function is calculated from the 
nucleon density p(b) distributed according to the Woods-Saxon formula (Equation (|2.17|l ): 

T(b) = J p(b,z)dz (2.22) 

Figure |2~. 101 compares the hard-sphere cut-off (at b = 2Ra u ) and P no hadronic- 

The probability of mutual nuclear excitation with breakup PxnXn was first com- 
puted by Baltz, Klein and Nystrand in [7j for the case of ultra-peripheral photo-nuclear p° 
production accompanied by a mutual nuclear dissociation. The computation assumes the 
flux of virtual photons to be distributed according the Weizsacker- Williams approximation 
(|2.8|) . The lowest-order probability for an excitation of a colliding beam ion to any state 
which emits one or more neutrons (Xn) is: 

[■ d 2 n (u, b) a, i 

Pkxn)(P) = / ) 9 -> (o,) a - (2.23) 

J d 2 b u 

where cr 7j 4-»A* (w) is cross-section of the excitation with breakup of nucleus A by a single 

photon of energy u. This quantity has been determined by measurement at a wide range 

of energies [5]. 

9 We distinguish between the case where both excited nuclei emit exactly one neutron (In, In) and the 
case where both excited nuclei emit one or more neutrons (Xn,Xn). 
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Figure 2.10: Probability of no hadronic interaction taking place between the ions as a 
function of nuclear impact parameter. Solid line - theta- function cut-off, dotted line - 
Pno hadronic (b) computed with Woods-Saxon nuclear form- factor. 



At small impact parameters, Pqixu)^ can excee d 1 and can be interpreted as 
a mean number of excitations. The unitarization procedure, similar to the unitarization 
for e + e~ pair production probability, needs to be utilized. The probability of having at 
least one Coulomb excitation is then Pc(Xn)(b) = 1 — exp ^— P^ X n) fii) ■ I n the reaction 
involving the dissociation of both ions, each individual breakup occurs independently |14j . 
The probability is thus the square of the individual breakup probabilities, i.e. PxnXn(b) = 

{P C (Xn)(b)) 2 . 

Figure 12.111 shows the probability of mutual Coulomb nuclear excitation with 
breakup computed for 120 values of b between 5 fm and 60 fm for Au ions at y / s/\F/v = 200 
GeV. For the excitation and breakup of Au ions to be likely, the impact parameter must be 
small (b < 30 fm). 



2.4 Other Related Processes at RHIC 

2.4.1 Ultra-Peripheral Electromagnetic Processes 

The bound-free pair production is another pure QED process in ultra-peripheral 
relativistic heavy-ion collisions which is of practical importance in the collider. It is the 
process, where a pair is produced but with the electron not as a free particle, but into an 
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Figure 2.11: Probability of mutual Coulomb nuclear excitation with breakup as a function 
of the nuclear impact parameter. 

atomic bound state of one of the nuclei. Using the Equivalent Photon Approximation and 
the approximate wave functions for bound state and continuum, the cross-section for this 
process for Au beams at RHIC was found to be about 90 bjllj. This changes the charge of 
the ions in the beam, causing beam loss. 

Electron and positron can also form a bound state, positronium. For the two- 
photon processes, the positronium can only be in the para-positronium state 1 Sq. An 
interesting state with a negative charge parity - orthopositronium ( 3 5i state of positronium) 
- can be created in three-photon interactions |41j. like the one in Figure l2~o1 Cross-section for 
this process is suppressed with respect to para-positronium production only by {Za) 2 and 
in the case of Au ions at RHIC the production of orthopositronium should be comparable 
to the production of para-positronium. 

A wide range of final states not involving e + e~ pairs is possible in the two-photon 
interactions at RHIC. The photons couple to any state with internal charge constituents 
(i.e. quarks) and spin/parity J PC = h ,0 ++ ,2 h or 2 ++ (for J < 2). The spin-1 final 
state in a two-photon reaction is impossible, because spins of massless photons can only be 
aligned or antialigned. We discuss some of the final states below, grouping them by the 
area of physics interest. These topics represent excellent future two-photon physics research 
opportunities at RHIC [EH E2 H2| . 
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• High- field QED: t + t and pairs 

These events are created by the same mechanism as e + e~ pairs. The cross-sections 
are low, due to the very high masses of these leptons. 

• Meson Spectroscopy: Search for glueballs and other exotica 

These events consist of a single meson or an exotic particle (e.g. glueball) in the 
final state. Because two photons couple to (charge) 4 , the cross-section (1(77 — ► X) 
is a direct measurement of the quark and isospin content of the mesonic final state. 
Consequently, final states consisting of charged particles ( qq ) are possible, but pure 
gluon final states are not possible. 

• Meson form factors (pairs) 

Two-photon interactions also produce meson pairs via 77 — » XX. At the hadron 
level, photons couple only to charged mesons, so it +/ k~ should be produced, but not 
7r°7r°. At high enough energies photons couple directly to the quark content of the 
mesons, and both tt + tt~ and 7r°7r° final states are produced in comparable numbers. 
By comparing the rates of the two final states, the transition can be studied, and the 
size of the mesons determined. 

2.4.2 Ultra-Peripheral Photo-nuclear Processes 

An ultra-peripheral collision of heavy ions might involve only one photon, interact- 
ing with the hadronic field of another ion, in which case a reaction is called 'photo-nuclear' 
interaction. One of the possible descriptions of this process is that a photon from one ion 
can fluctuate into a virtual qq pair, which scatters on the other ion and emerges as a real 
vector meson. 10 For Au ions at RHIC design luminosity this reaction produces several vec- 
tor meson species: p° at the rate of ~ 100 Hz, u at the rate of 12 Hz and <j> at the rate of 
7.9 Hz |53| . 

The p° is a short-lived resonance (width T = 150 MeV) which quickly decays into 
a 7r + 7r _ pair. Additionally, 7r + 7r~ can be created in RHIC collisions via the non-resonant 
channel 7 — ► 7r + 7r~ . These processes have been recently observed at RHIC by the STAR 
collaboration at -^/saF/v = 130 GeV and y / s/v r /v = 200 GeV both as an exclusive channel and 

10 This interaction mechanism is called 'vector dominance model', and believed to be dominant at RHIC 
energies. There are a number of other approaches to photo-nuclear processes, including photon-Pomeron 
scattering. 
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Figure 2.12: STAR collaboration observed ultra-peripheral n + Tr~ pairs. Left: pair trans- 
verse momentum distribution. Right: da/dM nn distribution (dashed lines represent direct 
7r + 7r~ and p° contributions) |24j . 

in coincidence with mutual nuclear excitations |24l H"5j . This is a potential background for 
e + e~ production, as discussed in Chapter d 

Figure 0. 121 shows the pair transverse momentum and invariant mass distribution 
of identified ultra-peripheral coherent events Au + Au — ► Au* + Au* + ir + ir~ at ^/snn = 130 
GeV. The transverse momentum distribution shows a prominent peak at low momenta. The 
incoherent tt + tt~ pair spectrum is modelled with the same-sign pair combinations (ir + ir + , 
ir~ir~, shaded histograms in Figure I2.12J) . and doesn't show a peak at low transverse mo- 
menta. The invariant mass distribution is split into contributions from the coherently 
produced p° mesons (Breit-Wigner distribution), direct 7r + 7r~ (fiat in Mi nv ) plus the inter- 
ference between these two sources and background from un- identified e + e~ pairs at very 
low invariant masses. 
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Chapter 3 



Monte Carlo Event Generator 



The following chapter describes a method we used for generating e + e~ with the 
nuclear excitation. The method is based on the Equivalent Photon Approximation. The 
output of the event generator will be used in Monte Carlo studies of detector acceptance, 
efficiencies and systematic effects (Chapter EJ). The events are generated in the STAR 
laboratory frame of reference. Figure 13.11 shows the STAR coordinate system definition, 
used throughout the rest of this thesis. 



3.1 Differential Cross-section Computation 

To generate e + e~ events, we need first to compute the differential cross-section 
da e + e - /dMi nv dY . As we see from the equation (|2.7j) this quantity depends on the two- 
photon luminosity with nuclear excitation. We perform numerical integration over two 
impact parameters b± and 62 in (x, y) plane to find the two-photon density (see Figure l3~Tj) : 
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(3.1) 

Equation 1)3. 1|) is derived from the original Equation l|2.9|) with the addition of the 
requirement of the simultaneous nuclear excitation (PxnXn) and replacement of a ^-function 
cut-off with a probability of no hadronic interaction between the ions (P no hadr)- The upper 
limits of integration are taken as the distance from the nuclei where the field becomes very 
small: i? max ~ A>y Au /u]. 
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Figure 3.1: Left: STAR coordinate system with a schematic view of the STAR TPC and 
the RHIC beams. Right: two-dimensional photon density integration in STAR coordinate 
system. 

The value of the two-photon density ()3.1|) is computed for 100 values of the pair 
invariant mass (Mi nv = 2^fuj\uj<i) and pair rapidity (Y = 1/2 log (wi/u^)) in the range of 
100 MeV < Mi nv < 300 MeV and \Y\ < 1.3. We chose these values because they represent 
the kinematical range in which the STAR detector has non-zero acceptance to the e + e~ 
pairs. All mass and rapidity values are equidistant, with the distances AMi nv = 2.0 MeV 
and Ay = 0.026. Each value in this two-dimensional table is then converted into the two- 
photon luminosity (Equation Q2.1UJI ). and multiplied by the lowest-order QED differential 
cross-section for the two-photon annihilation into a e + e~ pair, taken at the corresponding 
values of the pair mass and rapidity (Equation (|2.6| ). As the result we get a tabulated 
differential distribution da e + e - / dMi nv dY . 

We expect that the use of EPA for cross-section studies of the e + e~ pairs with 
high invariant masses is justified, since the photon virtuality is significantly smaller than 
the e + e- pair invariant mass (Q 2 ~ l/R 2 Au ~ (30) 2 MeV 2 < M 2 nv ). The use of Rau as a 
lower intergration limit in (|3.1|) is also justified by this circumstance. 

3.2 Event Generator 

We generate electron/positron tracks (track momenta p e + and p e -) in the de- 
tector lab frame. The origin of the tracks is taken to be at (0, 0, z), where z is generated as 
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a Gaussian random variable with mean zero and a = 30 cm to approximate the distribution 
of vertex position in the data. 1 The following describes how we draw p e + and ~p e -- 

3.2.1 Drawing Y and M inv 

Once we have obtained a tabulated differential cross-section distribution we can 
draw two correlated random variables Mi nv and Y according to this distribution. A two- 
step procedure needs to be used, whereby the value of Mi nv in the range (M^ 1 , M^^) = 
(100 MeV, 300 MeV) is drawn according to its marginal distribution, and then the value of Y 
in the range (y min , y max ) = (0, 1.3) is determined according to the conditional distribution 
f Y (Y'\M inv ). 

We draw a random number from a uniform (0, 1) distribution and obtain a value 
of comparing this random number to a table of values of 

TUT ymax 

/ / dYdM'^j^ 

^^min ymin 

(3.2) 



i nv 
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da 
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inv 

Then this value M® nv is used to determine Y° by comparing another draw from 
uniform (0, 1) distribution to a table of values of 

Y 

J U1 dM inv dY 

(3-3) 

1 max 

f dY — 

J al dMi nv dY 

^min 

and randomly assigning the sign of Y . 2 Using these values of Y° and Mf nv the momenta 
of the interacting photons can be calculated. 

3.2.2 Drawing Transverse Momentum of the Photons 

We draw a value of the photon transverse momentum after we determine the 
photon energy (for EPA photons \p z \ = E). Since the transverse momentum is much 

lr The smearing of the z vert in the data is significant (a z w 30 cm, see Chapter[lJ. Therefore the simulations 
include z vert smearing to study its possible effects on the track reconstruction. The smearing of the vertex 
radial position is small {or ~ 0.25 mm) and is expected to have negligible effect on the reconstruction. We 
do not include transverse vertex smearing in the simulations. 

2 We use the property that the distribution of Y is symmetric around zero. 
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smaller than the energy of the photon, this is a reasonable approach. Adding the transverse 
momentum after the photon energy has been determined makes photons slightly virtual, 

2 2 9 *? 9 

since (q^) = (E,p x ,p y ,p z ) = E — p^_ — E < 0. For a given value of photon energy a> 
the photon transverse momenta pj_ are distributed according to the distribution (|2.19|) . It 
is impractical to use the same method as we used for M{ nv and Y to draw a value of p± 
since a new table of the distribution function F p± (u) would need to be generated for each 
new value of uj. Instead, we use a Von Neumann accept/reject method |64| to make a single 
random draw from a distribution F p± (u) for a photon with energy u>. 

The x and y projections of the photon momentum are obtained using a randomly 
chosen (uniformly distributed) azimuthal angle. The transverse momentum of the two- 
photon system p 1 ^ is the sum in quadratures of the transverse momenta of the two individual 
photons. 

At this stage we can find the momenta of the photons in their center of mass 
system, by boosting the lab photon momenta by a Lorentz-factor 

~P={p t ° t ,p t y ot ,p t z ot )/(E 1 + E 2 ) (3.4) 

Due to non-zero transverse momenta, the two-photon rest frame is not collinear 
with the z cm -axis. The Euclidean rotation which yields the frame collinear with photon 
momenta is specified by two parameters: the axis of rotation ~u and the angle of rotation 



tot ^ 3 

l J cm A A cm 5 



) (3-5) 

3.2.3 Generating Angular Distribution of the e + e Pairs 

We start by generating electron and positron track momenta in the e + e~ center 
of mass frame. The angular distribution for the pair (in the frame of reference collinear 
with photon momenta) is given by (|2.5jl . We use accept/reject method to draw the value 
of 9*, and we draw azimuthal angle of the pair ip from a uniform (0, 2tt) distribution. To 
transform the momenta into the frame of reference which is collinear with the lab frame, we 
z ~z = (0, 0, 1) in the CM frame. 
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Generated Total Pair Transverse Momentum, GeV 




Total Transverse Momentum, GeV 

Figure 3.2: Total pair distribution for pairs with 100 MeV < M inv < 300 MeV and 
\Y\ < 1.3. 

perform a Euclidean rotation of the generated momenta around axis ~u by the angle —<f>, 
defined in Equation (|3.5j) . 

The final step is to transform the electron and positron track momenta from the 
center of mass frame into the lab frame. This requires a Lorentz boost by — f3 (defined in 
Equation (|3.4j0 . Appendix El provides formulae for the boost and rotation transformations. 

3.3 Monte Carlo Results for AuAu -> Au*Au* + e+e 

Kinematical Distributions 

A few interesting observations can be made at this point. First, we now can plot 
the distribution of the total transverse momenta of all ultra-peripheral e + e~ pairs. This 
spectrum is a convolution of the photon energy spectrum and the energy-dependent single 
photon transverse momentum distribution. Figure 13.21 shows this distribution for AuAu 
collisions at ^/snn = 200 GeV. The distribution displays a prominent peak at p 1 ^ ~ 5 
MeV/c, which is a defining signature of ultra-peripheral reactions at heavy-ion colliders. 

We also present the distribution of the generated pair invariant mass and the pair 
rapidity (Figure l3.3|) . The invariant mass spectrum falls off very quickly with increasing 
Mi nv . Total pair rapidity distribution is shown with a theoretical prediction (dashed line) 
between -3.2 and 3.2. The distribution reaches maximum at zero, and drops off for high 
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Generated Pair Invariant Mass, GeV Generated Pair Rapidity 




Pair Invariant Mass, GeV Rapidity 



Figure 3.3: Left: generated invariant mass spectrum. Left: pair rapidity spectrum, dotted 
- theoretical prediction, solid - generated. 



Electron Rapidity vs. Positron Rapidity, |Y|<0.1 
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Electron Rapidity vs. Positron Rapidity, 0.9<|Y|<1.1 
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Figure 3.4: Contour plots of electron vs. positron rapidity for pairs with \Y\ < 0.1 (left) 
and 0.9 < \Y\ < 1.1 (right). 



values of \Y\. 

Unlike the total pair rapidity, the individual track rapidities are centered away from 
zero rapidity, as shown in Figure l3~H For a pair with a very central rapidity (|Y| < 0.1) 
the individual track rapidity reaches a maximum at ~ 1.5 and reaches up to 5.0. Typically, 
electron and positron tracks have individual rapidities of similar absolute value, but of 
opposite signs. For the pairs with non-central rapidity (0.9 < \Y\ < 1.1) the track rapidity 
peaks at ~ 2.0 and reaches up to 6.0. Since both the electron and positron tracks are 
ultra-relativistic, the individual track rapidities (Y e + and Y e -) are extremely close to the 
individual track pseudorapidities (r/ e + and r/ e -). 

Figure 1531 shows a dependence of the total pair rapidity Y distribution on the cut 
on the maximal absolute value of individual track pseudorapidity. The cuts [|^ e +| < r/ max , |^ e -| < r/ max ] 
effectively limit the total pair radity: \Y\ < ^ max . The suppression is strongest for the larger 
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values of |y|. 



Pair Rapidity Distribution, | r| track | <1.0 
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Figure 3.5: Effect of the r] max cut on the pair rapidity Y distribution. Solid: generated Y 
distribution (scaled), dashed: Y distribution after ry max cut 
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Track p x : M lnv > 0.1 GeV, | n eleclron | <1.15, 1 r, po8llton |<1.15 



Track Pl : M,„„ >0.15 GeV, | n eteelr0 „ |<1.15, | n po , illon |<1.15 
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Figure 3.6: Single track transverse momenta for pairs with Mi nv > 100 MeV (left) and for 
pairs with Mi nv > 150 MeV (right). Track pseudorapidities are limited to \>q track \ < 1.15. 



It is also important to examine the individual track transverse momentum spectra 
for tracks at mid-rapidity (\rj e \ < 1.15). Figure l3~^l shows the spectra for such tracks with a 
cut applied to the lowest invariant mass of the e + e~ pair. We see that applying a cut on 
the invariant mass limits from below the individual track transverse momentum. 



Approximating 9* by 9' 

Since the event generator retains the information about the generated angle 9* for 
each event, we can compare these values to the values of angle 9' computed solely from the 
observed momenta of e + e~ tracks. Figure l3~Tl compares the two distributions. We see that 
the variable 9' approximates 9* very well in the range of 45° < 9* < 135°, and there is a 
significant difference between the two variables for 9* ~ 0° and 9* ~ 180°. This is due to 
the fact that the photon polar angle (9 p hoton ~ VtAu) cannot be neglected in comparison 
with 9* if 9* « 0° or 9* « 180°. 

Histogram of 9 * (solid) and 9 ' (dashed) Histogram of ( 6 • ■ 8 ') , 45 ° < 9 * < 135 




120 140 160 180 

Angle, deg 




Figure 3.7: Comparing 9* vs. 9' for generated e + e pairs. Left: ID histograms of 9* (solid) 
and 9' (dashed). Right: 9* - 9' for 45° < 9* < 135°. 
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Cross-section Predictions in the Limited Kinematical Range 

The Monte Carlo routine we described can provide a prediction of the e + e~ pair 
production cross-section within a limited kinematical range. This is important since the 
STAR detector has a limited acceptance for individual charged tracks in p t ^ ack and rf rack . 
Table l3~T1 shows several Monte Carlo cross-section predictions for e + e~ production at RHIC 
with nuclear excitation in a limited kinematical range. 



Cuts 


a (mb) 


10 MeV < M inv 


5- 10 4 


100 MeV < M inv 


195 


100 MeV <M inv and \Y\ <3.0 


142 


140 MeV <M inv <265 MeV and |Y|<1.15 


46.40 


140 MeV <M inv <265 MeV, |Y|<1.15 and cos(6>')<0.7 


4.95 


140 MeV <M inv <265 MeV, |F|<1.15 and \ v track \ < 1.15 


2.58 


140 MeV <M inv <265 MeV, |F|<1.15 , \ v track \ < 1.15 and p*_[ acfc > 65 MeV/c 


2.08 



Table 3.1: e + e production with mutual nuclear excitation cross-sections for various 
acceptances. 



35 



Chapter 4 

Experimental Apparatus 

RHIC is a new accelerator facility in Brookhaven National Laboratory which 
started functioning in the summer of 2000. It can collide head on protons and a wide 
range of heavy nuclei at very high energies. 

There are four heavy-ion experiments at RHIC named BRAHMS, PHOBOS, PHENIX 
and STAR placed at four beam intersection regions at the collider. Each of these detectors 
is optimized for measuring different final states, but there are overlaps in their capabilities 
so that consistency checks can be made between them. STAR detector has the advantage 
of the complete 2tt azimuthal angle coverage over the central rapidity region. Each of the 
RHIC intersection regions is also equipped with two identical Zero Degree Calorimeters for 
monitoring the collisions. 

4.1 RHIC System 

The RHIC accelerator complex reflects a long history of collider / accelerator devel- 
opment at Brookhaven National Laboratory. Figure POl shows the layout of the accelerator 
facility, starting from the production of a gold beam at the Tandem Van de Graaf to accel- 
eration to full energy at the main RHIC ring. 

At the beginning, partially ionized gold atoms are emitted from a source, such 
as a high temperature gold filament. The positively charged ions are accelerated through 
the Tandem Van de Graaff's two 15 million volt electrostatic accelerators, and are passed 
through thin sheets of gold foil, which further ionize the gold atoms. Ions exiting the 
tandem enter the Alternating Gradient Synchrotron (AGS), where they are accelerated in a 
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257-meter diameter radio-frequency synchrotron to a total energy of 11 GeV/nucleon. The 
AGS employs focusing technique with its 240 magnets situated along the acceleration ring. 
Each magnet successively alternates its magnetic field gradient inward and outward from 
the ring, allowing the beam to be focussed in both horizontally and vertically. The final 
complete ionization of the gold ions also happens at the AGS. Prom the AGS the ions are 
diverted to a transfer line to the main RHIC ring, where a switching magnet injects ions 
into a counterclockwise and clockwise rings. 

The main accelerating rings at RHIC known as a yellow (clockwise) and blue 
(counterclockwise) rings, are each about ~ 610 m in radius 25j. The rings are filled by the 
ions in a boxcar fashion, resulting in 57 distinct bunches of ions (each containing ~ 10 9 
197 Au +79 ions) in each ring. The acceleration of bunches in these synchrotron rings is 
achieved by radio-frequency cavities and bending the beam into a circular shape by super- 
conducting magnets, positioned around the ring. Using the super-conducting technology 
RHIC is able to collide the heavy ions at the center-of-mass energies higher than any other 
machine in world (table l4~Tj) . The top center-of-mass energy achievable at RHIC for heavy 
ions is ^/snn ~ 200 GeV per nucleon in the lab frame. 

Once the desired collision energy is reached the beams are synchronized to cross 
in six interaction regions. The design frequency at which RHIC bunches crossed, which 
depends on the energy of the beams, was 9.37 MHz in 2001. In the interaction regions, the 
beams are focussed and steered by quadrupole magnets for collisions at approximately 180 
degree angle (head-on collisions). The design length of the region of space where collisions 
take place ('interaction diamond') is about 20 cm. In the 2001 run the interaction diamond 
had a (almost) Gaussian profile, with a sigma of ~ 60cm. At top energy, 200 GeV, in 2001 
RHIC has achieved AuAu luminosities of ~ 10 26 cm~ 2 s —1 |30j . 

The length of time that RHIC can continuously provide collisions at a single time 
(called 'fill') is limited by a few factors. First is the beam loss. This is mostly attributed to 
the collisions of the beam ions with the residual atoms inside the vacuum beampipe, called 
the beam-gas collisions. The next most significant source of beam loss was the capture by the 
positive Au ions of an electron from an e + e~ pair, discussed in Section[23I Another process 
which contributed to the beam loss, though to the lesser extent, was single photonuclear 
excitation of the beam ions followed by nuclear breakup [23} . The synchrotron radiation 
was not a significant source of beam loss, due to the very high masses of the accelerated 
ions. The second reason for interrupting the fill is the slow dispersion of the bunches in the 
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Year 


Facility 


Location 


Species 


Efoeam 


Ecm 


1974 


- 1991 


BEVALAC 


LBNL 


197 . 197 . 

Au+ Au 


2 GeV 


2 GeV 


1994 


- present 


AGS 


BNL 


197 . 197 . 

Au+ Au 


11 GeV 


5 GeV 


1994 


- present 


SPS 


CERN 


208 „. 208 „, 

Pb+ Pb 


158 GeV 


17 GeV 


2000 


- present 


RHIC 


BNL 


197 . 197 . 

Au+ Au 


100 GeV 


200 GeV 



Table 4.1: Some beam species and collision energies for various heavy-ion accelerators. The 
facilities are given by a facility acronym and location. Beam energies correspond to the top 
per nucleon beam energy achievable by the facility. The Ecm column shows the per nucleon 
energy of the system in the center-of-mass frame jS2j . 

longitudinal to the beam direction. This causes the interaction diamond to become very 
long, and a large fraction of collisions unusable. The typical fill length was four to six hours 
in the 2001 run. 

Relativistic Nuclear Physics Program at RHIC 

The central focus of the RHIC physics program is study of nuclear matter at high 
temperatures and densities. In the central collision of the highly relativistic nuclei enormous 
energy densities (> lGeV/fm 3 ) are created. These conditions may create a system of the 
theoretically proposed state of deconfined quarks and gluons - a Quark Gluon Plasma [Glj. 1 
This system expands rapidly and the temperature and density will drop below critical 
values causing the formation of the gas of interacting hadrons. As the hadronic system 
expands and dilutes, inelastic interactions between hadrons become more and more scarce. 
At the point these inelastic interactions stop, the system reaches a so-called chemical freeze- 
out [S3]. Eventually, the system expands so much that even the elastic interactions between 
the hadrons cease, and the system is said to have reached a thermal freeze-out. RHIC was 
specifically designed with the goal to observe the novel state of matter QGP and to study 
the exact details of freeze-out of this state into a normal hadronic matter. 

Another fundamental question addressed by RHIC is a spin structure of nucleon. 

The total spin- 1/2 of a proton is the sum of the spin contributions of the quark constituents 

of the hadron, their angular moment and the gluons. In the present understanding of 

the nuclear spin, quarks contribute only 1/3 of the total nuclear spin, and the gluon spin 

1 Om Universe may have been in the state of the Quark Gluon Plasma a few microseconds after the Big 
Bang. 




Figure 4.1: RHIC system. STAR is at the point labeled 2 [2H]. 
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contribution is non-negligible 18., but no measurement exists yet. RHIC is planning to 
collide polarized beams of protons to study the details of the gluon spin term. 
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Figure 4.2: Left: Location of the ZDCs around the interaction point. Right: hardware 
components of the ZDCs |23| , 

Zero Degree Calorimeters 

The main devices for monitoring the rate of RHIC collisions are Zero Degree 
Calorimeters. These detectors are located on either side of each of the six interaction 
regions at RHIC, 18 meters from the interaction regions along the beam pipe. In STAR, 
the two ZDCs are labeled 'East ZDC and 'West ZDC The ZDCs are capable of detecting 
and measuring the energy of uncharged particles emitted from the interaction region in the 
direction of the beam [23;. The detectors consist of layers of Cherenkov fibers sandwiched 
between the tungsten absorber plates (10 x 10 cm in size in the direction transverse to the 
beam, tilted at 45° angle with respect to the beam), as shown in Figure [Ql When hadrons 
hit the detector hadronic showers are developed in tungsten absorbers, and generate a signal 
in the Cherenkov fibers. 2 The signal from the fibers is then sent to the photomultiplying 
tubes (PMTs), with the summed analog output of the PMTs forming the ZDC signal (ZDC 
ADC signal). 

The detectors have a 98% ± 2% efficiency (flat over all detector area) for detecting 

uncharged particles within a 2 mrad cone around the beam direction [23]. This acceptance 

region is sufficient to detect all of the spectator neutrons in AuAu collisions. Such neutrons 

2 This is a slight modification of the traditional sampling hadronic calorimeter design, which uses scintil- 
lators for hadronic shower sampling. 
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are typically moving at the speed of the beam in the direction of the beam, and their 
transverse momentum p± is mostly due to the Fermi motion (~ 5 MeV/c)jH]- 3 Thus the 
angle between the direction of the beam and the momentum of the neutrons is of the order 
of 1.4 mrad, which is within the acceptance. 

The ZDCs are read out every RHIC bunch crossing and provide information 
whether the crossing resulted in a collision or not. The detectors also provide the timing 
information, measuring the arrival time of the neutral products of the heavy-ion collision 
to the ZDC on either side of the interaction point. In effect, the arrival time difference in 
two ZDCs provides a measurement of the primary interaction location in the direction of 
the beam. The resolution of the method can reach a = 3.2 cm|23j. 

4.2 STAR Detector 

The Solenoidal Tracker at RHIC (STAR) shown in Figure fOl is a detector designed 
primarily for measurement of hadron production over a large solid angle. STAR features a 
detector system for high precision tracking, momentum analysis, and particle identification 
at the mid-rapidity. The large acceptance of STAR makes it particularly well suited for 
event-by-event characterization of heavy-ion collisions and for reconstruction of particle 
decays, such as K s — ► tttt and A — > pir 32 . 

The detectors in STAR can be roughly split into two groups - fast detectors, which 
can read out data at the frequencies close to the frequency of RHIC bunch crossings, and 
slow detectors, which operate significantly under the RHIC frequency, but can provide a 
much more detailed information. A discussion of each of the detectors used in the 2001 run 
follows. 

4.2.1 STAR Magnet 

The main STAR detector magnet is a 1100 ton solenoidal structure that consists 

of a 130-turn aluminum solenoid conductor. The magnet structure encloses most of the 

other STAR detectors, as shown in Figure [Ql It can provide 0.5 T of magnetic field when 

energized to a current of 4500 A, and 0.25 T with the current of 2250 A|26j. The field flux 

is returned through the poletips at either end of the solenoid, and then through a set of 

3 For Mutual Coulomb Dissociation neutrons p± distribution is even more narrow than Fermi 
distribution 23 . 
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Magnet 






Figure 4.3: Schematic view of STAR detector [3*2*] . 



flux return bars, which are situated on the cylindrical surface of the solenoid. The purpose 
of the flux return is to keep the field uniform within the volume of the solenoid. 

The STAR magnet serves two purposes for the detection of charged products of 
the heavy-ion collisions at RHIC. First, it allows the momenta and sign of charged particles 
to be calculated by measuring the curvature of the particle track as it passes through the 
field volume. Second, the B-field produced by the STAR magnet is oriented in the beam 
direction, as is the drift field used in the TPC. As we explain in the following section, this 
reduces the dispersion in ionization produced in the TPC, increasing the position resolution 
of the TPC. 

4.2.2 STAR Time Projection Chamber 

STAR Time Projection Chamber (TPC) is a main tracking detector in STAR. It is 
a large gas filled chamber capable of measuring three-dimensional space-points along charged 
particle trajectories. In comparison to silicon detectors, the TPC has coarser position 
resolution, but can make multiple measurements over a large volume. TPC is an intrinsically 
slow detector, since its readout speed is always limited by the drift speed of the ionized gas 
atoms ( ~ 40ug)pHj. 
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Figure 4.4: Cylindrical geometry of the STAR TPC's gas volume and the sector layout on 
the endcaps. The boundaries of the TPC gas volume are the Inner Field Cage (IFC), Outer 
Field Cage (OFC) and the TPC endcaps [SSj. 

TPC Main Gas Volume 

The STAR TPC gas chamber has a cylindrical geometry which extends 4.2 m in 
length and 2 m in radius, as shown in Figure [Ql The ionization region or active volume of 
the TPC is more than 45 m 3 . This volume is kept slightly above atmospheric pressure (2 
mbar) and filled with 10% CH4 and 90% Ar gas (P10). Signals originate from electrons that 
are freed when moving charged particles ionize the gas. The positive ions and free electrons 
move apart under the influence of a 147 V/cm electric field between the central membrane 
and end caps of the TPC. The positive ions are carried to the cathode at the central 
membrane and the electron clouds drift towards the ends of the detector. Positive ions are 
neutralized when they reach the cathode plane and the electron clouds are amplified in a 
Multi-wire Proportional Chambers (MWPC) close to the end caps. Since the drift velocity 
of the electrons is known, one coordinate (z) of the starting point can be deduced from the 
time taken for the electrons to drift to the MWPC. The other two coordinates are found 
through the projection of the signal onto a pad plane mounted below the MWPC. The 
pad plane lies perpendicular to the beam axis and is segmented into 136,608 pads. The 
electronics are capable of recording 512 time bins from each pad, of these, about 348 are 
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Figure 4.5: The three wire planes of the Multi-wire Proportional Chambers of the STAR 
TPC. The outer and inner sector geometries are shown on the left and right, respectively 

E5|- 



read out between the central membrane and the MWPC. In total, the volume is effectively 
divided up into more than 47 million space-points. 

The Multi-wire Chambers shown in Figure 14.51 have three planes of wires. The 
anode wires are closest to the pad plane and are 20 fim in width. The inner sector anode 
wires are set to 1170 V and the outer to 1390 V. The combination of the fine width and 
high voltage on the anode wires produce a strong radial field near the surface of the wires. 
Drift electrons create ionization avalanches as they accelerate towards the positive anode 
wires. Positive ions created in these avalanches produce image charges on the pad plane. 

The ground plane grid is the middle plane of wires that separates the drift volume 
from the amplification region. The ground plane grid has three main functions: to provide 
a ground plane for the drift field, to shield the pad plane from the gated grid and to capture 
some of the positive ions created near the anode wires. The drift field is established between 
the -31 kV central membrane and ground grid plane. The ground plane grid significantly 
reduces the signals induced on the pads when the gated grid opens. This prevents these 
induced signals from compromising the resolution of ionization signals at the beginning of 
the drift period. A large fraction of the slowly drifting positive ions created near the anode 
wires are neutralized on the wires of the ground plane grid. Positive ions that drift into the 
active volume, leakage current, cause distortions in the drift field. 

The gated grid is furthest from the pad plane. The main purpose of the gated 
grid is to stop non-triggered ionization from reaching the amplification region and stop 
positive ions created in the amplification region that leak past the ground plane grid from 
reaching the active volume. In the closed state, adjacent gated grid wires alternate from 
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Outer Pads Inner Pads 

6.2 mmx 19.5 mm 2.35 mmx 11.5 mm 




Figure 4.6: Pad plane of one TPC sector. Each sector contains 5692 pads. The inner and 
outer pad geometries differ to compensate for the radially decreasing hit density |35| . 

positive to negative potentials. These potential differences set up electric fields between 
the wires that are perpendicular to the drift direction. The fields capture the non-triggered 
electrons and positive ions. In the opened state, the voltage on the gated grid wires is set 
to the corresponding equipotential surface of the drift field. In this state the gated grid is 
transparent to the drift electrons. 

The TPC pads are laid out in sectors that cover 30° in azimuth, as shown in I4.4I 
There are 24 identical sectors mounted on the east and west ends of the TPC. Each sector 
has 13 inner and 32 outer pad rows, as shown in Figure l-Ofl Effectively, the pads are plate 
capacitors. Local electric field changes are created on the surface of the pads by the slowly 
drifting positive ions created in avalanches near the anode wires. These local field changes 
induce currents on the pads and subsequently in the TPC electronics. 

TPC Electronics and Data Acquisition 

The TPC readout electronics boards are mounted on the back of each sector. Each 
sector has 181 analog Front End Electronics boards (FEE) and six digital readout (RDO) 
boardsE36j. The circuitry on each FEE is separated into two parallel 16 channel circuits and 
is capable of covering up to 32 pads. The analog signals on the TPC pads are amplified, 
shaped, stored and digitized in two chips on the FEE. The Pre-Amp/Shaper-Amp (SAS) 
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amplifies and shapes the signal. The SAS feeds the 512 slot switched capacitor array (SCA). 
This chip is an analog storage unit that also contains an analog to digital converter. The 
chip allows fast, low-noise sampling of the signal with minimal power consumption. It 
also permits digitization and readout of data at a reduced rate. Upon request the SCA 
chip digitizes the stored voltages on the capacitors and passes them onto a multiplexer 
on the RDO board. The multi-plexer communicates via fiber optic links with the data 
acquisition (DAQ) crates $1 . DAQ Receiver Boards (RB) receive the data from RDO 
boards and perform pedestal subtraction, zero-suppression and charge cluster finding along 
the TPC time dimension. The data is then sent into the DAQ Event Builder (EVB), which 
is responsible for gathering the event fragments from all of the sectors of all detectors, 
packing the data into DAQ files, buffering the data on disk and then sending the data to 
the High Performance Storage System (HPSS). 

4.2.3 Additional Detectors 

A number of other detectors were present in STAR in the 2001 run, but not used 
in this analysis. Here is a brief description of them. 

Forward Time Projection Chambers fFTPC) |33|. located on both sides of the TPC, 
provide track momentum reconstruction in the forward region (pseudorapidity interval of 
2.5 <| 7] |<4.0). These detectors are especially useful for asymmetric collision studies, such 
as p+A collisions, and any other events where tracking in the forward region is desired. 
The 2001 run was a commissioning run for these detectors. 

The Central Trigger Barrel fCTB) j4"0] is a cylindrical detector positioned just 
outside of the outer diameter of the TPC, as shown in Figure 1131 The slats of the barrel 
are connected to photomultiplier tubes which give a response proportional to the number 
of charged particles which interact with the scintillator medium inside the slats. Thus CTB 
is a tool for measuring the charged multiplicity in the central rapidity region of | ij |< 1.0. 
The CTB can be read out every RHIC bunch crossing (every 104 ns). 

The Barrel Electromagnetic Calorimeter (EMC) [37! is also located just outside 
of the outer field cage of the TPC. The detector is a barrel made up of 120 modules of 
a lead scintillator sampling calorimeter. This detector is intended to study rare, high-pj_ 
processes and photons, electrons, 7r° and rj mesons in the same rapidity region as the CTB. 
The detector can be read out on the nanosecond scale. In the 2001 run the EMC was 
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commissioning the first set of modules. 

The Multi-wire Proportional Chambers (MWPC) 56 is a detector which is a part 
of the STAR TPC, that can also be used as an independent charged multiplicity counter in 
the region of 1.0 <| rj |< 2.0. This detector is located on the endcaps of the STAR TPC, and 
can provide a pixelized count of the charged particles traversing its active volume for each 
RHIC strobe. The MWPC was expected to be of great use for studying ultra-peripheral 
events with forward-peaked tracks, and a STAR Note 434:56 describes a simulations pack- 
age to estimate its efficiency for detecting such tracks. Due to noise, the MWPC could not 
be used as an independent detector in 2001. 



The Silicon Vertex Tracker (SVT)jHHl is a detector positioned right around the 



interaction region, outside of the beam pipe and inside the Inner Field Cage of the TPC. Its 
main goal is to provide additional tracking in the area directly adjacent to the interaction 
region, thus greatly increasing STAR tracking resolution and vertex finding accuracy. 



The Ring Imaging Cherenkov Detector (RICHEST] was placed just outside the 



CTB. The RICH covers a small area, 1 m x 1 m and is designed to provide high precision 
velocity measurements for the high momentum particles that pass through it. This enhances 
the particle identification capabilities of STAR for high momentum tracks. 



The amount of material in the path of a charge particle can be expressed in terms 



of a radiation length[BI|. 4 Table IO summarizes the thickness of material (expressed in 
radiation lengths) traversed by the track between the interaction region and the active 
volume of the main tracking detector TPC. 



The traversed radiation length defines the amount of secondary scattering suffered 



by a particle. The total angle of the secondary scattering over the length of the track can 
be approximated as a Gaussian random variable with zero mean and the width: 



4.2.4 Material Table in STAR 




(4.1) 



4 Radiation length is a mean distance over which a high-energy electron loses all but 1/e fraction of its 



energy. 

5 For events with vertex position within \z vert \ < 75 cm. 

6 For a particle traversing a straight line between 50cm and 100cm in radial direction. 
7 Does not affect the scattering of detected particles inside the TPC. 
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Structure 


Radiation Length (%) 


RHIC Beampipe 5 


0.28 


SVT 


6.00 


IFC 


0.62 


TPC gas 6 


0.39 


OFC 7 


1.26 


Total without OFC 


7.29 



Table 4.2: Material traversed by a particle in STAR. Table compiled from [SHj and |39j . 

where p is a particle momentum, (3 - velocity, Z is the charge of the particle and Xq is the 
traversed radiation length. Deviations from the particle's trajectory as it travels through the 
material in STAR can be calculated using Equation (|4.1[) and for a 100 MeV/c particle will 
be on the order of 0q ~ 0.37 rad. This will cause mis-reconstruction of the scattered track 
momentum and limit the precision with which the track momenta can be reconstructed by 
STAR. 

4.3 Triggering in STAR 

As we mentioned before, STAR can read out and store on tape the data at the 
rate of 100 Hz, where as the RHIC bunch crossing frequency (maximal possible collision 
rate) was 9.37 MHz. This is a frequency at which the RHIC RF cavities operate, and this 
frequency is supplied to all experiments at RHIC, called a 'RHIC strobe'. We need a set of 
detectors which are capable of reading out information at the frequency of the RHIC strobe 
and making a rough decision if a particular bunch crossing might contain an interesting 
collision. Such detectors are the CTB, the EMC, the MWPC and the ZDCs. In the 2001 
run only the CTB and the ZDCs were used for triggering. 

4.3.1 STAR Trigger System Design 

The on-line analysis of trigger detector information (or 'triggering') in STAR is 
done in several steps, called Levels through 3 of the trigger. Each increasing level of trigger 
has access to a more detailed information from the trigger detectors and correspondingly 
takes more time. The whole trigger system is synchronized with the RHIC strobe, so that 
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the trigger detectors are read out and the data is passed between the different levels of the 
trigger system only when the beams cross in the interaction region. Below we discuss the 
STAR trigger design in the 2001 run. 

Hardware Triggers (LO - L2) 

Level is the basic hardware trigger layer. This layer looks at every RHIC crossing, 
deciding whether to accept the event or not. Level consists of a multi-layer system of 
data storage and manipulation boards (DSMs) [2*7] . Each DSM consists mainly of a field 
programmable gate array (FPGA) which can be configured by using the VHD-Language. 
For each RHIC strobe, trigger detector data is passed into the first layer of the DSMs, from 
where the data is pipelined into the following layers of DSMs. The physical algorithms 
implemented in the FPGAs reduce the data along the DSM layer structure, yielding a 
final trigger decision, which combines the data from the trigger detectors (e.g. charged 
particle multiplicity count in CTB). In the final DSM layer this information is combined 
with detector LIVE/BUSY signal from the slow detectors and a trigger is issued if the slow 
detectors are live. 8 The time allowed between the RHIC crossing and Level decision is 
1.5 microseconds. 

Trigger Level 1 uses CPUs to analyze the output of the first layer of the DSM tree 
during the TPC drift time. Level 1 has 100 microseconds to either abort the event pass it 
to Level 2. Level 2 also consists of CPUs which have access to the full trigger information 
(first-layer DSM inputs, most finely grained information). The Level 2 analysis takes place 
during the TPC digitization time (~ 5 milliseconds). If Level 2 doesn't abort the event, the 
event is sent to the data acquisition system. 

On-line Event Reconstruction (Level 3 Trigger) 

The Level 3 trigger is a processor farm which can use the information from the 

STAR TPC to perform an on-line reconstruction of events at the rate of ~ 100Hz"2*5*j. Given 

a reconstructed event, Level 3 is able to analyze the properties of this event (such as total 

charge in the event, number of tracks in the STAR TPC, track dE/dx, etc) and make a 

decision if this event should be written out on tape or not. Thus Level 3 trigger can select 

events based on physics observables, such as rare particles like J /\P or antinuclei. This 

8 'Live' is a term which means that the slow detectors are finished processing the previous event and are 
ready to accept a new event. 
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trigger is particularly useful for selecting events with low multiplicity and specific event 
topological signature, such as ultra-peripheral collisions. 

4.3.2 Types of Triggers Available in STAR 2001 Run 

Given the trigger capabilities described above, a few different trigger conditions 
were programmed into the STAR trigger logic, called 'trigger types'. We discuss the trigger 
types that were used for studying ultra-peripheral collisions. 

The Minimum Bias trigger was programmed at Level to require a coincidence 
signal in the East and West ZDCs (more than ~ 40% of a single neutron energy deposition in 
both ZDCs). The dominant fraction of events collected with this type of trigger was hadronic 
collisions with charged multiplicities in STAR TPC between ~ and 3000, since nearly all 
hadronic interactions of the colliding beams result in the emission of spectator neutrons in 
both directions of the beam. Additionally, a fraction of ultra-peripheral collisions can occur 
simultaneously with the photonuclear excitation of both Au ions, resulting in a Minimum 
Bias trigger. 

The next modification of the trigger was to accept at Level events which satisfy 
the Minimum Bias trigger conditions and show a collision that has a vertex position in 
the direction of the beam (determined by the ZDCs) within 30 cm from the center of the 
interaction region. This trigger type was called a Minimum Bias Vertex trigger. The reason 
for this type of trigger is that we want the events that happen in the center of the STAR 
TPC, where the geometrical acceptance is symmetrical in the direction of the beam. 

Another trigger type was called a 'topology trigger' and was designed specifically 
to trigger on 2-track ultra-peripheral events. This trigger was programmed in the Level of 
the STAR trigger system to accept events which had charged multiplicity 1 in both north 
and south quadrants of the STAR CTB, but nothing in the top and bottom quadrants (to 
ensure cosmic ray events would not be accepted). Optionally, this trigger could also utilize 
the capacities of Level 3, requiring that events have zero total charge and have a vertex 
position in the beam direction within 100 cm of the center of the interaction region. 

An important feature of STAR trigger system is that it can simultaneously look 
for events satisfying either one of several different trigger types. This is called parallel trig- 
gering. This allows a more complete use of the RHIC run time, since several different kinds 
of physics events can be collected in the same run. For instance, for a significant portion 
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of the 2001 run the Minimum Bias and topology triggers were run in parallel, allowing the 
STAR physicists to collect both ultra-peripheral data and hadronic collisions |27j. 
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Chapter 5 

Event Reconstruction and 
Simulations 

The following Sections describe the off-line event reconstruction in the STAR TPC, 
with a special emphasis on low-momentum electron tracks. With a few exceptions, this 
reconstruction chain also applies to the simulated events, which we discuss in the Section 

5.1 Reconstruction 

In order to extract meaningful physics from the raw data, the pixel information 
which is stored from each tracking detector must be 'converted' into reconstructed tracks. 
This is implemented with event reconstruction software, organized into a chain of recon- 
struction procedures (reconstruction 'makers'). In this analysis the only tracking detector 
used is the TPC, and the main purpose of the reconstruction software is to determine charge 
clusters in pad time-space and convert them into position coordinates in the TPC. Using 
the cluster information, pattern recognition software finds tracks which are in turn used 
to accurately determine the position of the beam interaction point (primary vertex). The 
tracks that are found at this stage are either 'secondary' (not coming from the primary 
vertex, such as cosmic muon tracks, for instance) or 'primary' (produced at the interaction 
vertex), so a re-fit of the tracks is undertaken, including the interaction point in the fit. 
Given the curvature of the track and the magnetic field, we can measure the rigidity of 
the particle - momentum divided by the charge. To get the momentum, tracking software 
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s>0 




Figure 5.1: Geometric definitions of some commonly used track parameters. The crosses 
mark possible data points. ^ is a azimuthal angle of the track momentum at the point 
(xo,yo). pt, p z and p are the track transverse, longitudinal and full momenta. A, R and 
are the helix dip angle, radius of curvature and phase at the point (xo>2/o) |6U| . 



assumes an absolute value of the charge equal to the charge of the electron. 

5.1.1 TP C Hit and Global Track Finding 

Electron-cluster finding along the TPC time dimension (z) is performed in DAQ 
crates, before the data is shipped to HPSS. The remaining hit reconstruction is accomplished 
off-line by scanning along rows for adjacent pad signals. Each pad row is treated separately 
and hits are never reconstructed across different pad rows. Thus, a track crossing the entire 
active gas volume of the TPC from the inner field cage to the outer field cage can have up to 
45 hits. A cluster is then a set of signal sequences on adjacent pads on a pad-row which have 
some overlap in time. The distributions are fit with a two dimensional Gaussian function to 
determine the centroid of the cluster along the pad-row and in the time direction. This is 
assumed to be a hit position. Once the 3-dimensional positions of all hits are determined, 
the hits are used by a TPC tracking algorithm, where a hit represents a possible point on 
a track. 

Tracking starts at the outermost pad-row of the TPC, where the hit density is 
lowest. It searches for groups of three hits that lie close in physical space. These groups 
are the initial "seeds" of tracks. Track segments are formed from seeds through linear 
extrapolations to adjacent pad rows and clusters in z. Since the STAR magnetic field 
guides charged tracks into helical trajectories, segments are extended inward and outward 
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with helical extrapolations. After these extensions, a minimum number of 5 hits is required 
on each track so that the five parameters of the track's helix model are uniquely defined. 
If this requirement is not met, then the track is flagged as a bad fit. Tracks are fit in two 
independent projections of the helix. A circular projection is fit in the plane perpendicular 
to the magnetic field, the xy-plane. A line is fit in the sz-plane, where s is the length along 
a track's circumference [00]. Tracks found with such a procedure are called 'global tracks'. 
Kinematic variables of a track are calculated from the parameterization of the helix; found 
via the circular and linear fits, and the direction and magnitude of the magnetic field. The 
definitions of some track parameters are given in Figure 15.11 and a more complete list of 
commonly used track parameters and their definitions can be found in Appendix ^ Each 
global track is saved in the computer memory along with its kinematical and geometrical 
parameters and the basic information about the track fit, such as the number of hits used 
in the track fit and the \ 2 °f the helix fit divided by the number of hits used. 

The tracks which make a full 360° rotation (in transverse projection) in STAR 
TPC will be reconstructed as two separate tracks, each making only a half-turn. A side- 
on view of such tracks is given in Figure 15.31 each of the two real tracks in the event is 
reconstructed as two half-spirals. One of the tracks will point towards the true interaction 
vertex in coordinate z (tracks 1 and 3), while the other will not (tracks 2 and 4). Tracks 
like the tracks 2 and 4 are called spiraling tracks. Spiraling tracks cause the contamination 
of clean 2-particle events with low-momentum tracks, as we will discuss in Chapter El 

5.1.2 Low Multiplicity Vertex Finder 

An accurate knowledge of the interaction position (primary vertex) for each col- 
lision is required as this helps distinguish between track which originate from the primary 
vertex (primary tracks) and those from backgrounds (cosmics, beam-gas interactions, etc) 
or from weak decays. The primary vertex coordinates can also be included in a track re-fit 
for primary tracks which leads to a better determination of their momentum. In this section 
we describe a low multiplicity vertexer (LMV) - a routine which is used to find primary 
vertices for events with low multiplicity. 

The vertex finder procedure starts with an assumption that the true vertex is 
located at (0, 0) in xy-plane. 1 Any track which does not come to this point closer than 

lr This is a well justified assumption, since typically the vertex position is only within a = 0.25 mm, from 
(0, 0) in the a:y-plane, and the ^-position of the vertex is within a — 60 cm from the center of the interaction 
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a certain minimal value Rmax is removed from the vertex finding routine. The remaining 
tracks are modeled as helices and extrapolated inside the IFC, and an analysis is run to 
determine the amount of multiple scattering suffered by these tracks in the TPC gas, the 
TPC IFC, the SVT and the beampipe. The more secondary scattering the track experienced, 
the less weight this track will be given in the vertex fitting. Finally, a primary vertex is 
determined as a point (xo,yo,zo) which minimizes the weighted sum of squared distances 
from each helix to this point. A check procedure is performed to find the helix which is the 
farthest from the found vertex, and to remove the track if that distance is greater than a 
pre-set maximum DCA max . If a track was removed from an analysis at this stage, a new 
vertex fit is performed, if all tracks were retained, the vertex finding routine is finished. The 
error squared of the vertex position is set equal to the weighted sum of squared distanced 
from the tracks to the primary vertex. 

There are a number of reasons why a vertex finding routine can fail for an event 
consisting of several low-momentum tracks which come from a true primary vertex within 
an interaction diamond. Fist of all, a track's kinematic parameters may be sightly mis- 
reconstructed, and the track may be found to pass within more than R m ax from the beam 
in the xy-plane. Secondly, after the vertex fit is performed, some tracks may be found 
to pass within more than DCA max from the primary vertex, this happens again due to 
mis-reconstruction of tracks' kinematic parameters. Since each of these occurrences results 
in a removal of an outlier track from the analysis, the number of tracks remaining for the 
determination of the vertex may be less than two, which automatically causes the vertex 
finding routine to fail. Thus, there is a certain probability that for the low multiplicity 
events consisting of low-momentum tracks a vertex finder will find no vertex, even when in 
fact there is a vertex. We will estimate this probability in Section f6. 2.41 

5.1.3 Primary Track Fitting and Kalman Filter 

Including the primary vertex in the track fits can significantly improve the accuracy 
of the track reconstruction. A global track pointing back to the primary vertex and having 
an extrapolated distance of closest approach to the vertex less than 3 cm is considered 
to have originated from the vertex. Global tracks which satisfy this criteria are termed 
primary tracks and are subjected to a re-fit of the track. The primary vertex is included 
diamond. 
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in the fit. Since the initial estimate of the tracks' kinematic parameters are know from the 
helix fit of the global tracks, a more sophisticated procedure, called Kalman filter can be 
utilized (see |46| for for a general description of Kalman filter). This is more realistic than 
a simple helix fit as it takes into account multiple scattering and energy loss in the inner 
field cage and the gas inside the TPC in each step of the fitting routine, in order to obtain 
the best estimate of the track momentum. This approach requires three passes through the 
data. The first pass, known as filtering, starts at the outermost pad-row and works its way 
in, removing hits from the track that fail very crude cuts. The second pass, referred to as 
smoothing, starts at the innermost radius and works its way outwards. Again, outliers are 
removed from the track but this time the rejection cuts are much harsher. The third pass 
is for the evaluation of the track parameters and the total x 2 - 

In order to estimate the degree of multiple scattering and energy loss, an estimate of 
the mass of the particle must be used. This is accomplished by assuming that all particles are 
pions. This is true for approximately 70-80% of tracks in a typical STAR high-multiplicity 
event, but obviously not correct for the current analysis. Although it is desirable to fit every 
track with a correct mass hypothesis, this was not performed in the 2001 run reconstruction, 
as the process is extremely CPU intensive. We will study the effect of the incorrect mass 
hypothesis on the reconstructed primary track momenta in Section \7. II 

5.1.4 Particle Identification with the TPC by Specific Energy Loss 

Charged particles lose energy while traversing matter. The main process contribut- 
ing to this energy loss for heavy particles is Coulomb scattering of the traveling particles 
with electrons in atomic orbits, causing ionization of the traversed medium. The signal 
heights of the pads in a hit cluster give a measure of the energy loss in the TPC gas. These 
signals are summed over all pads in a cluster. The energy loss per unit length for each hit is 
then equal to this sum divided by the length of the track over the pad row. The distribution 
of the energy loss samples for hits along a track has a Landau-like distribution 64] . An es- 
timate of the mean of this distribution is calculated by computing the average of the lower 
70% of the distribution - this procedure, called truncation, is necessary since it ensures 
that the truncated mean is not sensitive to the extremely large values of dE/dx in a small 
fraction of clusters. The result is a measure of the energy loss per unit distance along a 
track, dE/dx (KeV/cm) in the TPC gas. Figure E2] shows the scatter plots of the dE/dx 
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Symbols 


Definition 


Units or Value 


m e 


electron mass 


0.51 MeV/c 2 


r e 


classical electron radius 


2.81 fm 


N A 


Avogadro's number 


6.022 xl(T 23 mor 1 


K 


AitN A r e 2 c 2 


0.307 MeVg^cm 2 


q 


charge of the particle 




z 


Atomic number of the medium 




A 


Atomic mass of the medium 




I 


Mean excitation energy 


MeV 


T 

max 


maximum kinetic energy that can be imparted to a free electron in a collision 


5 


Bethe-Bloch density effect correction factor 



Table 5.1: Physical constants used in the Bethe-Bloch function, Equation (|5.1|) . 



distribution vs. the momentum of the tracks, measured by STAR TPC. 

For moderately relativistic charged particles with m 3> m e the ionization energy 
loss is a dominant process, and can be described by the Bethe-Bloch function pi]: 

- (dE/dx) truncated = Kq -— I -In j 2 (3 - - \ (5.1) 

The variables used in the Equation (|5.1|) are defined in table 15.11 According to 
Equation (|5.1|) . ionization loss only depends on a particle's charge (q), velocity {(3) and the 
properties of the medium that is traversed. To get the expected mean dE/dx as a function 
of the particle momentum in STAR TPC we need to assume a mass and a charge for a 
particle. Figure IS~21 shows in solid lines the resulting functions for several particle mass and 
charge assumptions. 

For lighter particles, such as electrons, several other processes are important - 
Bremsstrahlung and Moller scattering (Bhabha scattering for positrons) are the two most 
significant. Bremsstrahlung is particularly important, since electrons are very light, and in 
the region of momentum of the order of dozens of MeV/c they are ultra-relativistic. The 
total length of the STAR TPC in the transverse direction (150 cm between the inner field 
cage and the outer field cage) is equal approximately to 1.17% of the radiation length. Thus, 
an electron traversing a TPC with a momentum of ~100 MeV/c will lose ~ 0.0117- (^-)- 100 



dE/dx Vs. P | 




Figure 5.2: Scatter plot of track dE/dx measured by the STAR TPC vs. track momentum 

m- 

MeV~ 0.7 MeV due to Bremstrahlung, which is comparable to the ionization energy loss 
over 150 cm of traversed P10 gas. The energy loss by Bremsstrahlung will (mostly) not 
be detected by the STAR TPC, therefore the Bethe-Bloch formula is not adequate for the 
description of energy loss by electrons in the STAR TPC. In general, we observe that the 
description of the energy loss by electrons in the STAR TPC is not well described in the 
GEANT detector simulation package (next section). 

Despite the aforementioned difficulties, for the electrons in the region of momenta 
from ~10 MeV/c to ~100 MeV/c, star TPC detects dE/dx energy loss which is clearly 
different from other species of particles (the lowest band on the left in Figure I5.2|) . The 
average dE/dx for the electrons shown as a solid line was not computed with the Bethe- 
Bloch formula (|5.1|) , but rather determined experimentally from a clean sample of electrons 
and positrons from photon conversions p2j. The electron band is sharply crossed by other 
particle bands as they approach their minimum ionizing values, as illustrated in Figure [5~2l 

With a high level of confidence particles can be identified in regions where the 
dE/dx bands are not overlapping. In regions where bands are partially merged, simultane- 
ous Gaussian fits to the bands in momentum slices are used to extract the probability that 
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a track is of a given particle species. With the tracking field at 0.25 T, the STAR TPC 
has a one sigma dE/dx resolution of about 8.2% of the dE/dx value. The knowledge of 
both a predicted dE/dx value and the resolution of dE/dx is a tool for performing particle 
identification. Using the predicted dE/dx, a measured dE/dx and a sigma, a known frac- 
tion of a certain particle band can be sacrificed in order to eliminate other particles. This 
is accomplished by requiring the dE/dx of all particles to fall within a certain number of 
sigma of a predicted band [dE/dx deviant). 

5.2 Details of the Simulation Chain 

We used the events generated with a Monte Carlo generator, described in the 
Chapter Next we used the GEANT detector description package 1 to simulate the 
passage of electrons and positrons from the collision point through the detector material 
and ionization of the TPC active gas volume by the charged tracks. GEANT also simulates 
the creation of secondary particles in the detector material interactions. 

The secondary particles consisted mostly of knock-out electrons and a small num- 
ber of e + e~ pairs produced by hard Bremsstrahlung photons emitted by the original elec- 
trons or positrons. These particles are added to the table of Monte Carlo generated tracks 
with unique ID's, separating them from the original event generator tracks (primary parti- 
cles) . 

The ionization produced by charged tracks traversing the TPC gas was passed 
through a TPC Response Simulator (TRS). TRS simulates the TPC response from the 
drift of the ionization in the TPC gas to the output of the front-end electronics. It includes 
the drift, diffusion, amplification and response of the electronics for the electrons created 
by ionization in the TPC gas. Simulated pad signals can then be then passed to the event 
reconstruction routine. 

Additional Procedures Applied to the Monte Carlo 

An additional step in the event reconstruction routine for the simulated events 
is the association of the reconstructed tracks with the Monte Carlo tracks in the STAR 
TPC. The association routine takes the TPC hits used for a fit of a given track and checks 
if those hits might come from any single Monte Carlo track. If more than 5 TPC track 
hits are associated to the Monte Carlo track, the reconstructed track is associated to the 
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Figure 5.3: Side-on view of the STAR TPC. Shown are two particles going through the 
TPC, each particle will be reconstructed as two separate tracks. 

Monte Carlo track. In the conditions of the low-multiplicity environment association has 
an efficiency of ~ 99% ±1%, nearly all reconstructed tracks (with number of hits per track 
> 5) were matched to the Monte Carlo generated tracks. We use only the reconstructed 
Monte Carlo tracks that are associated to the primary particles in Monte Carlo efficiency 
studies. 

Also, we wish to exclude from Monte Carlo efficiency studies the spiraling tracks. 
Figure IB~3l shows two particles in the STAR TPC, one in a positive z direction, one in the 
negative. Each one of the tracks is reconstructed as two separate tracks. To exclude the 
second parts of the spirals (tracks 2 and 4) we apply the following requirement: 



where Zt rac k is the z position of the first innermost hit on the track (shown as z\ , z2 , z3 
and z4 in the Figure IS~3|) . z ver t ex is the z position of the GEANT generated vertex, and 
Vtrack is the reconstructed pseudorapidity of the track. Monte Carlo studies show that the 
requirement (|5,2j) eliminates 99% of the spiraling tracks. Aditionally, we determined that 
for simulated e + e~ pairs with a reconstructed primary vertex, spiraling tracks make up 
only a 2% contamination. 

5.3 Event Pre-Processing and Data Storage 

After the event reconstruction, all of the data is written to a Data Summary Tape 
(DST). The DSTs are split into chunks of approximately 130 events as each file is large 
(~ 500 megabytes) by this stage. Each event is a separate entity on the DST, with a list of 




(5.2) 
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all primary and global tracks in the event, a map between the primary and global tracks, 
a vertex information, and general characteristics of the event, such as the time the event 
was taken, the trigger type for this event, the conditions of the STAR detector and so 
on. Additionally, if the event comes from the simulation, the full information about the 
simulated event kinematics and the association map between the simulated tracks and their 
reconstructed counterparts is included. 

This specific analysis utilized one additional stage of data processing, unique to 
the ultra-peripheral events. Such events can typically be easily separated from the hadronic 
events by either the low multiplicity of tracks in the event, or if an event was triggered 
with a specific ultra-peripheral type of trigger (such as the topology trigger, Section 14.31 
for instance). All DST files were filtered for the ultra-peripheral events based on the above 
conditions, and all the selected events were analyzed by a specific ultra-peripheral analysis 
routine. The routine takes all tracks and pairs them into all possible pair combinations. For 
each pair several characteristic variables are computed, such as the total charge of the pair, 
the total transverse and longitudinal momentum, the invariant mass of the pair assuming 
both tracks are electrons, pions or kaons, and a few more. This routine was run separately 
on the primary and global tracks in each event. The pairs comprised of primary tracks were 
called 'primary pairs', and the pairs comprised of the global tracks were called 'secondary 
pairs'. 
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Chapter 6 

Analysis 



In order to compute the cross-section of a reaction, we need to count the number 
of events which are believed to come from a studied reaction (iV b S erved events); know the 
total integrated luminosity of the data set considered (Ltot) and the efficiency of triggering 
on the events of this type (Eft* 1 ^), and to determine the efficiency for reconstructing the 
events of this type (Eff dctector ) : 

-^observed events /g <\ 

a ~ L tot x Efr* rig x Eff dotector 
This chapter describes the selection cuts we used for finding events of the type 
AuAu — ► Au* Au* + e + e~ in the pool of reconstructed events in the STAR detector, and 
the computation of total luminosity times the trigger efficiency in the analysis dataset. 
The subsequent chapters (0 |HJ) will discuss the reconstruction efficiency and the measured 
cross-section. 



6.1 Data Set 

We desire to observe the ultra-peripheral electromagnetic production of e + e~ pairs 
with simultaneous nuclear excitation in gold-gold collisions: AuAu — > Au*Au* + e + e~ . 
Events of this type constitute a sub-set of the events collected with Minimum Bias Trigger 
at 0.25 T magnetic field (half- field) . During the 2001 data taking the largest half-field 
dataset was collected with a Minimum Bias Vertex trigger, this is the data set we used 
for the analysis. The dataset has about 800,000 events, about 60,000 of these events were 
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found to have less than 10 global tracks per events. The presented analysis will consider 
these events. 

Luminosity 

The Equation 16.11 can be applied equally to ultra-peripheral reactions and to 
hadronic reactions. In particular, this equation can be applied to a total hadronic cross- 
section for gold-on-gold reactions, which has been calculated to be 7.2 b ±10% [2]. Simul- 
taneous photonuclear excitations of both Au ions contribute additional 3.2 b cross-section 
in the Minimum Bias dataset 23 . Analysis of the hadronic reactions in STAR detectors 
shows that events with more than 7 negatively charged hadrons in the STAR TPC in the 
pseudorapidity range \rj\ < 0.5 and with transverse momenta p± > 100 MeV/c constitute 
80% of the total hadronic cross-section [30]. This allows us to measure the total luminosity 
times the trigger efficiency in the Minimum Bias Vertex data sample ('normalization to 
hadronic cross-section'): 



L . Eff trig - — MinBiasVertex dataset (8+ positive hadrons in \rj\ < 0.5 and pj_ > 100 MeV/c) 
tot ~ 0.8 x 7.2 b 

(6.2) 

In the Minimum Bias Vertex data set, Ltot ■ Eft*" 8 for ultra-peripheral events can 
be determined by normalization to the hadronic cross-section only for events which have a 
longitudinal vertex position within \z ver t\ < 25 cm. This is due to the fact that to use nor- 
malization to hadronic cross-section, we need to make sure that the trigger efficiency is the 
same for the hadronic events (which are usually high-multiplicity and produce a lot of neu- 
trons in the ZDCs) and ultra-peripheral collisions (which are low multiplicity and typically 
produce just a few neutrons in the ZDCs). As we mention in Section [4.11 the ZDCs have 
close to 100% efficiency for detecting both single and multiple neutrons; therefore Minimum 
Bias trigger efficiency is the same for hadronic and ultra-peripheral events. However, the 
Minimum Bias Vertex trigger uses ZDC timing information for rejecting a portion of events 
that are accepted by the Minimum Bias trigger, and the efficiency of this rejection may 
be different for the events with a lot and a few neutrons. Selecting only the events within 
\ z vert\ < 25 cm, which is 5 cm tighter than the ZDC timing vertex cut, 1 allows us to avoid 
the problem of unequal trigger efficiencies for hadronic and ultra-peripheral events. 
Vertex position resolution in ZDC timing cut is on the order of 3.2 cm ("Section 14. II . 
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The luminosity times trigger efficiency for the data set available to us with a cut 
\ z vert\ < 25 cm was found to be 62 mb _1 , with systematic uncertainty of 10% due to 
uncertainties in total hadronic cross-section. 

6.2 Event Selection 

We will be working with the sub-set of the Minimum Bias Vertex triggers, a total of 
about 60,000 events (with less than 10 global tracks per event). We will begin by identifying 
criteria for selecting ultra-peripheral events among the other low-multiplicity events. We 
then present a Monte Carlo study of the STAR detector acceptance for low-momentum 
tracks, and restrict our attention to the tracks which are in the region where STAR detector 
has good acceptance. Selecting e + e~ pairs among other ultra-peripheral events requires 
particle identification with dE/dx. Finally, we present a study of vertex finding efficiency 
and restrict the position of the primary vertex for selected events. 

6.2.1 Identifying Ultra-Peripheral Pairs 

We require that the total pair charge be zero. If this requirement is not met the 
pair is guaranteed to be a background, and we will take advantage of this property later 
(Section 17. 3j) . We also require that there be exactly one primary pair in an event, and 
at most two more non-primary tracks - this allows for additional tracks due to occasional 
contamination by beam-gas interactions, spiraling tracks and secondary particles. 

A crucial way of distinguishing the pairs produced in ultra-peripheral coherent 
interactions of two gold ions is to cut on the maximal value of the total transverse momentum 
of the pair. The physics motivation for this cut was explained in Sections 13.31 and 13.2.21 
Figure 021 shows that the total transverse momentum of the coherent e + e~ pairs peaks at 
about 5 MeV/c and drops off significantly at 70 MeV/c. We set the (j/^) max cut at 100 
MeV/c for the coherent pairs. 

6.2.2 Track Acceptance Cuts 

The STAR TPC has a limited acceptance for charged tracks. First of all, in 
transverse plane tracks must have a diameter greater than the radius of the Inner Field 
Cage (Rifc = 50 cm), which translates into a requirement on the tracks' transverse mo- 
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Figure 6.1: Left: TPC tracking efficiency as a function of track p± for tracks within the 
\r]\ < 1.15 cut. Right: acceptance as a function of track 77, for tracks with p± > 0.65 GeV/c. 



menta p± >p m i n = 0.3qBRjpc = 37.5 MeV/c. Secondly, since the TPC has a finite length 
[Ltpc = 4 m), tracks must have pseudorapidity |r/| < — In (tan (0.5 tan" 1 [Rjfc/^-^Ltpc])) 
2 to make it inside the TPC. 

To find the regions of transverse momenta and pseudorapidity where the TPC 
tracking efficiency is relatively flat we used the Monte Carlo event sample described in 
Chapter |21 (individual track transverse momenta in the range 50 MeV/c < p± < 200 MeV/c 
and pseudorapidity \r]\ < 4 ). The tracking efficiency can be defined as: 



^ ^single tracks reconstructed, associated to primary particles ^ 

"' ' 111 ' ^single primary particle tracks generated 

The resulting distributions in single track p± and 77 are shown in Figure 16.11 

Figure l6~T1 demonstrates that STAR TPC has close to 1 acceptance for tracks with 

high transverse momentum at mid-rapidity. A slight dip in reconstruction efficiency as a 

function of rj at zero is due to the fact that tracks with \rj\ < 0.3 are more likely to have low 

p± and not be reconstructed than tracks with 0.3 < \rj\ < 0.7. A check procedure applying a 

cut pj_ > 0.95 GeV/c to Monte Carlo tracks, which ensures that reconstruction efficiency in 

p± is flat, returns an efficiency in rj which is flat for \rj\ < 0.7. We chose to place acceptance 

cuts on the primary track p± and r] at 

Ir/il < 1.15 M < 1-15 

(6.4) 

pi_L > 65MeV/c p 2 _L > 65MeV/c 
Figure shows that the TPC has some residual acceptance outside of the accep- 
tance region (|6.4[) . however we chose to cut away the regions of low acceptance, where the 
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Figure 6.2: Left: distribution of "Nhits in the data (dots) and Monte Carlo (solid black). 
Right: TPC tracking efficiency for tracks with N£™ = 6, N$£ = 8, N££ = 10, N££ = 12, 



simulation might not adequately describe the detector effect on the data. 

Figure 16.21 (left) shows the comparison of the number of hits per track (Nhits) 
distributions in the data with the cuts in Section T6 . 2 . 1 1 and Monte Carlo, and the agreement 
is quite good. We use Monte Carlo to study how the tracking efficiency depends on N hits' 
Figure 16.21 (right) shows the tracking efficiency vs. transverse momentum computed for 
tracks with more than 6,8,10,12 or 14 hits per track. For tracks with p± > 65 MeV/c, 
dependence of the tracking efficiency on the number of hits is negligible. STAR event 
reconstruction software writes on tape only TPC tracks with more than 10 hits per track, 
therefore we set N™ t ° = 10 cut in the analysis. 

6.2.3 Particle Identification via Specific Energy Loss 

Single track dE/dx depends on the track's momentum, particle mass and charge 
and in a lesser degree, on the number of hits on the track. To eliminate the dependence 
of track dE/dx on these variables, a new variable Z is introduced. The definition depends 
on a species of particles this variable is defined for. For pions Z n = log (^ ^E/dx) X (p) )' 2 ^ 
this formula is applied to a pion track dE/dx, the variable Z is independent of the track 
momentum and length. A set of data plots of track Z n for various slices of track momentum 
is shown in the Figure lo31 

In the Figure pions are on the right, electrons/positrons are on the left, except 
for the bottom right histogram, where pion peak is located on the left. Shape and the mean 
2 For a discussion of the expected value of dE/dx for electrons and pions see Section f5.1.4l 
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Z„(p) for track momenta .11 < p < .12 GeV/c 




Z„(p) for track momenta .1 2 < p < .1 3 GeV/c 




Figure 6.3: Z ff histograms, in p slices. Solid curves - Gaussian fits to the pion peaks. Tracks 
to the right of the dotted lines (upper histograms) are rejected as pions. 



of the pion band stay constant in all momentum slices. For a single track with momentum 
below 130 MeV/c, pions can be rejected at 95% confidence level by selecting tracks with 
Z n < —0.3. For the tracks with momenta between 130 MeV/c and 140 MeV/c rejection of 
pions at 95% confidence level cuts away a large fraction of electron/positron tracks, therefore 
we chose not to do particle identification for this p-slice tracks. For track momenta between 
140 MeV/c and 265 MeV/c electron and pion bands overlap completely. Thus, for tracks 
with momenta between 130 MeV/c and 265 MeV/c particle identification based just on the 
information from a single track dE/dx is impossible. However, doing a simultaneous e + e" 
dE/dx identification of e + and e~ tracks allows us to resolve this problem partially. For 
tracks with momenta above 265 MeV/c electron/positron tracks can again be separated 
from pion tracks by requiring > 0.25. 

The basic idea of the simultaneous pair identification is that if one track is identified 
as an electron (positron) with a high degree of likelihood, and the pair is identified as ultra- 
peripheral, the other track is most likely a positron (electron), since there are no processes 
that could produce a pair consisting of an electron and a charge-1 particle other than 
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a positron in the ultra-peripheral AuAu collision. 3 Table 16.11 explains the details of the 
simultaneous e + e~ dE/dx identification for a pair of tracks labeled track 1 (always a lower 
momentum track) and track 2 with momenta p\ and p2 respectively. 



track momentum configuration 


accept as e + e if 


1 


p x < 130 MeV/c, p 2 < 130 McV/c 


< -0.3and Z% < -0.3 


2 


pi < 130 McV/c, 13(Kp2 < 265 McV/c 


Z* < -0.3 and Zl within la of {dE/dx)™ ean 


3 


130<£>i<265 MeV/c, 130<p 2 <265 MeV/c 


identification via dE/dx impossible 


4 


pi < 130 McV/c, p 2 > 265 MeV/c 


Z*<-0.3 and Z%>0.25 


5 


130< pi < 265 MeV/c, p 2 > 265 MeV/c 


Zl > 0.25 and Z\ within 2a of (dE/dx)™ ean 


6 


pi > 265 McV/c, p 2 > 265 McV/c 


Z\ > 0.25 and Z^ > 0.25 



Table 6.1: Conditions that track 1 Z\ and track 2 Z^ must satisfy for a simultaneous dE/dx 
e + e~ pair identification. 



According to the table IB~T1 simultaneous e + e~ pair identification is possible when 
either one of the tracks has momentum below 130 MeV/c or above 265 MeV/c. As we dis- 
cussed in ChapterOl the number of produced e + e~ pairs drops very rapidly with increasing 
invariant mass of the pair (equivalently, with increasing track momenta). For this reason, 
the number of pairs with a track having a total momentum above 265 MeV/c is expected 
to be extremely small, and in the sample of 60,000 events preliminary analysis did not find 
any such e + e~ pairs. Therefore, to identify a pair as an e + e~ we require that the pair 
tracks' dE/dx satisfy either the requirement 1 or 2 in the table EU 

The requirement that at least one of the tracks should have a momentum less than 
130 MeV/c and the other track should have the momentum less than 250 MeV/c effectively 
limits the maximal invariant mass of the identified e + e~ pairs. We need to use Monte Carlo 
simulations to study the efficiency of simultaneous e + e~ identification as a function of pair 
mass. 

Pair Identification Efficiency from Monte Carlo 

We begin the pair identification study by comparing the electron track Z n dis- 
tributions in STAR TPC for the data and simulations. Figure IB~H shows the distribution 

Contribution from ultra-peripheral fJ> + fi~ pairs, where one muon decays via the channel fi — > e + vv is 
negligible due to the small fi + fi~ cross-section. 
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of track Z v for 3 track momentum slices. There is a significant difference in the mean 
and the width of the real electron peak (solid black histogram, left peak is electrons) and 
the simulated electrons (dashed histogram). This is due to the fact that GEANT package 
has an imprecise model for the energy loss by electrons. We applied mean shift and the 
width scaling to the simulations to make sure that the corrected simulated electron Z^ peak 
matched the data: 

. data I \ 
mean \ O ^ \P) 



zcorrected simu 



(p) = (zr u ( P ) - (zr u ( P )) mean ) ■ + {zt ata (?)) 



(6.5) 



where (Z% ata (p)) mean , ^ mu ( p )) mean , a f ta (p) and a s z imu (p) are the means and widths of 
Gaussian fits to the Z n data and simulations distributions in momentum slices 0.11 GeV/c 
<p< 0.12 GeV/c and 0.12 GeV/c < p < 0.13 GeV/c in Figure El 



Z„(p) tor track momenta 0.11 < p < 0.1 2 GeV/c 



Z„(p) tor track momenta 0.12 < p < 0.13 GeV/c 





UP) 



Figure 6.4: Comparison of data and GEANT Z n (p). Solid histogram - data, dashed his- 
togram - simulations. Data histograms show electron peaks (left) and pion peaks. Electron 
peaks in the data and simulations were fit with Gaussian profiles to find mean and width 
of the peaks. 

The pair identification efficiency can be computed as: 



pair identification 



^generated pairs passing cuts 1 or 2 



(6.6) 



^generated pairs 

The identification efficiency distribution as the function of pair mass and total 
pair transverse momentum is shown in Figure 16.51 The pair identification efficiency is 
close to 1.0 for low invariant masses (both tracks are in the region p < 130 MeV/c, where 
electrons/positrons are distinguishable from pions), and drops to zero for invariant masses 
on the order of ~ 300 MeV. We chose to place an upper cut-off on the pair invariant mass: 
Mi nv < 265 MeV to avoid the region where the efficiency is very low and the Monte Carlo 
does not provide an adequate description of the data. 
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Figure 6.5: Efficiency of pair identification vs. pair mass and total transverse momentum. 



Vertex Finding Efficiency 
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Figure 6.6: Vertex finding efficiency as function of pair invariant mass and transverse 
momentum. 



6.2.4 Vertex Finding 

Reconstruction of an e + e~ pair consists of track reconstruction of both tracks and 
finding a common vertex for these tracks. We need to study the vertex finding efficiency as 
a function of pair kinematical variables to ensure that we are not trying to detect pairs in 
the region where the vertex finding efficiency is low. 

We define vertex finding efficiency as: 



^reconstructed e + e pairs with a found vertex 



(6.7) 



^reconstructed e + e pairs within acceptance 
Figure I6.7I shows the vertex finding efficiency as a function of pair invariant mass 
and total transverse momenta. The efficiency is flat as a function of and Mi nv for 
Mi nv > 200 MeV, but drops off for invariant masses below 200 MeV. This is due to the 
fact that the low invariant mass pairs are composed of low momentum tracks, and the 
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Distance of Closest Approach between Global Tracks | | Radial Position of the Primary Vertex 




DCA, cm Radius, cm 



Figure 6.7: Distribution of the distance of closest approach between the global tracks and 
the vertex radial position. 

low momentum tracks undergo a lot of multiple scattering in the beam pipe or the inner 
field cage. As the result, the momentum of these tracks is often mis-reconstructed and 
projecting the tracks back to the common origin yields a very large distance of closest 
approach between the tracks, which causes the vertex finding routine to fail. We set a cut 
on the minimal invariant mass of the reconstructed pair at Mi nv > 140 MeV. 

We wish to restrict our attention to the events with the good quality of vertex 
finding. For the events with only two tracks used for the determination of the vertex, the 
quality of the vertex fit depends on the distance of closest approach between the global 
tracks which are used for the vertex finding {DCA). We limit this variable to less than 
4 cm. As a secondary check, we make sure that the found vertex is not too far (in the 
transverse plane) to the most likely position (0,0). This is the second cut: R ve rt < 9.0 
cm. The distribution of DCA and the vertex radius R ver t are shown in Figure loTTl for pairs 
which satisfy the definition of ultra-peripheral events (Section I6.2.1|) . To increase event 
sample, particle identification was not applied, but we selected events with both tracks' 
momenta below 200 MeV/c. This ensures that R ver t and DCA distributions are close the 
distributions for a pure e + e~ sample. The agreement between the Monte Carlo and the 
data is excellent, which means that the simulation of the multiple scattering in the detector 
and its effect on tracking and vertex finding is very accurate. 

The distributions of R ver t in the data and Monte Carlo have long tails, with ~ 3% 
of the evens having a reconstructed R ver t > 9 cm. This is due to the fact that determination 
of the point of the closest approach between 2 helical tracks of similar curvatures and 
opposite transverse momenta is unstable. For instance, in the limiting case of electron and 
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data (dots). Events with \z ver t\ < 25 cm 
Figure 6.8: xy projection of the STAR are used for luminosity normalization to 

TPC. Solid and dashed line arches repre- hadronic cross-section, 

sent electron/positron tracks. Points 1, 
2 and 3 represent possible reconstructed 
transverse vertex positions. 

positron tracks having the total transverse momentum of zero, the position of the point 
of closest approach between the two tracks in the transverse plane is completely undefined 
(Figure I6.8JI . However, the uncertainty in the R ver t determination has little effect on the 
determination of tracks' absolute momenta and the opening angle of the pair. 

Longitudinal Vertex Position Distribution 

Figure 16.91 compares the z ver t vertex distribution for ultra-peripheral events in 
the data and e + e~ Monte Carlo simulation. The Monte Carlo simulation (solid) has a 
Gaussian shape with half-width of 30 cm. The distribution of the vertices in the ultra- 
peripheral events in the data (dots) is close to Gaussian in the region | ^^er-t I < 25 cm, but 
drops of steeper than the Gaussian distribution outside of this cut. This is the effect of the 
vertex position cut in the Minimum Bias Vertex Trigger. In the data, the ratio of the events 
with a vertex inside < 100 cm to the events inside \z ver t\ < 25 cm is 1.52 ± 0.07. 

In the Monte Carlo, 99% of the events are within \z ver t\ < 100 cm. We have performed a 
check procedure on Monte Carlo simulation, grouping events in different slices of the z ver t 
vertex position, and the distribution of the kinematic variables of the reconstructed events 
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Wcis found, to be identical for all z ver t slices. 



Total Available Luminosity 

To increase the statistics, we would like to use ultra-peripheral events with a 
z vertex inside a 100 cm range, instead of 25 cm (see Section l6.1|) . For such events, the 
luminosity times trigger efficiency can be estimated from the luminosity for ultra-peripheral 
events satisfying a \z vert \ < 25 cm cut as: 

L tot ■ Eff^s = (L tot ■ Eff^s) . , _ x ^ l<10 ° (6.8) 

v ' \Zvert <25 AT. \ ' 

J " z|<25 

The total available luminosity was found to be 93.8 mb _1 with a systematic error of 10%. 
6.2.5 Summary of All Cuts 



Table l|6.2|) summarizes the analysis cuts. The first part of the table 16.21 contains 
cuts that select events with one primary e + e~ pair, produced in an ultra-peripheral elec- 
tromagnetic interaction. The second part of the table selects only those of the identified 
e + e~ primary pairs which are produces within the defined acceptance region. The third 
part selects the pairs which have high quality of the vertex finding. 



6.3 Resulting Raw Distributions 

After applying the cuts in table lfi~2l to the Minimum Bias Vertex dataset, we get 
52 identified events of the type AuAu — ► Au*Au* + e + e~. Figure IB. 101 presents an event 
display with x, y and z projections together with a side-on view of one of the identified 
events, showing the reconstructed tracks inside a STAR TPC. The event shows the typical 
characteristics of the reactions of this type: the tracks have a rather low transverse momenta 
(as evidenced by the strong curvature of the tracks) and are forward-peaked, with one track 
exiting the STAR TPC through the endcap. In the transverse projection (upper right figure) 
the tracks appear to be back-to-back and have very similar absolute values of the transverse 
momenta, this leads to a very low total transverse momentum of the pair, which is one of 
the distinctive characteristics of the e + e~ pairs produced in the coherent ultra-peripheral 
electromagnetic interactions of the heavy ions. 

Figure 16 . 1 II presents the raw (not corrected for detector efficiency) distributions of 
the pair kinematic parameters in the selected 52 events of the type AuAu — > Au* Au* +e + e~ . 
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Figure 6.10: Four views of an e + e~ pair in the STAR TPC. Top right: zy projection, top 
left: xy projection (transverse), bottom left: zx projection, bottom right: 3-dimensional 
view. 
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Identification Cuts 



Event triggered with Minimum Bias Vertex trigger 
Number of global tracks in event < 5 
Number of primary pairs in event = 1 
Primary pair total charge = 

Total transverse momentum of the primary pair < 100 MeV/c 



The z position of the primary vertex is within 100 cm from the zero 



EITHER pi < 130 MeV/c, p 2 < 130 MeV/c 


and Z\ < -0.3 , Zl < -0.3 


OR pi < 130 MeV/c, 130 <p 2 < 265 MeV/c 


and Zl<-0.3 , Zl within 2a of (dE/dx)~ 



Acceptance Cuts 
1 7/1 1 < 1.15 , rli > 65 MeV/c , \r] 2 \ < 1.15 , p± 2 > 65 MeV/c 
140 MeV < M inv < 265 MeV 



Data Quality Cuts 

Number of hits per track greater than 9 

Radial position of the vertex R ver t within 9.0 cm from (0, 0) 

Global counterparts of the primary tracks have DC A < 4 cm 



Table 6.2: Analysis cuts to select events of the type AuAu — > Au*Au* + e + e . 



The p^ distribution shows a clear peak at the low transverse momenta, which is a signature 
of ultra-peripheral e + e~ pairs. The distributions of <p and ip are flat, as expected for 
77 — > e + e~. The agreement between the Monte Carlo results and the data is good, which 
will allow us to use the Monte Carlo simulations for the efficiency corrections. 

Figure 16.121 shows the distribution of transverse momenta and pseudorapidities 
of electron tracks and positron tracks. The distributions agree well with Monte Carlo 
and within available statistics no difference can be observed between electron/positron 
kinematical distributions, suggesting that Coulomb corrections are not observed. 
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Figure 6.11: Raw kinematic distributions of identified e + e~ pairs in STAR, a) Mi nv , b) Y, 
c) p 1 "*, d) cos (9'), e) ip, f) Dots - data, solid lines - Monte Carlo. Dashed lines in items 
a) and c) represent acceptance region cuts. 
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Figure 6.12: Individual track p± spectrum (left) and r\ spectrum (right). Triangles 
positrons, crosses - electrons, solid line - lowest-order QED simulation. 



ZDC Spectra 

We do a check of the ADC signal spectra in the East and West ZDCs for the 
selected events. The signals in the ZDCs are from mutual electromagnetic excitations of 
Au ions emitting one, two or more neutrons. Figure 15. 131 shows the distribution of ADC in 
the East and West ZDCs and the fits of the data to the sum of two Gaussians, representing 
one neutron peak and two neutrons peak: 



f(x) = Aexp 



(x - x - sy 



(7- 



+ B exp 



(x — 2xo — sy 

2^2 



(6.9) 



To reduce the number of free parameters in the fit, only parameters A and B were found 
with a fit to the AuAu — > Au* Au* + e + e~data. The values of pedestal offset (5), first 
Gaussian mean (xo) and width (a) were taken from the analysis of ultra-peripheral p° 
production with mutual nuclear excitation 45 . 

Both histograms and their Gaussian fits show that 1 neutron emission events can 
be clearly separated from multiple neutron events by selecting events with ADC count 
below 27 in the East ZDC and below 25 in the West ZDC. Figure 16 . 141 shows the correlation 
of the East vs. West ZDC signals in the selected 52 events. 3 events satisfy the cut 
[ ZDC East < 27, ZDC West < 25] - these are (In, In) events. The ratio of (In, In) events 
to (Xn,Xn) events in the sample (0.06 ± 0.04) agrees within available statistics with the 
ratio of (In, In) events to (Xn, Xn) in the ultra-peripheral p° study (0.095±0. 05)115]. Since 
e + e~ production with nuclear excitation and p° production with excitation have the same 
mechanism of electromagnetic nuclear excitation, this confirms that the observed signals in 
the ZDCs are from mutual nuclear excitation decays. 
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East ZDC Energy Deposit (a.u.) West ZDC Energy Deposit (a.u.) 



Figure 6.13: ADC signal spectra in the ZDCs for events of the type AuAu — > Au*Au* + 
e + e~. Dots - data, solid lines - Gaussian fit. 
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Figure 6.14: Correlation of East ZDC ADC and West ZDC ADC for events of the type 
AuAu — > Au*Au* + e + e~ . Dashed lines separate 1 neutron peaks from multiple neutrons 
peaks. 
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Chapter 7 

Efficiency Corrections and 
Background Subtraction 

We use Monte Carlo simulations to do the efficiency correction for the observed 
spectra and to estimate the backgrounds. This chapter also describes the extrapolation 
technique we utilized to extrapolate the experimentally accessible cross-section to the 47r 
cross-section within the kinematical range of 140 MeV < Mi nv < 265 MeV and \Y\ < 1.15. 

7.1 Kinematical Parameter Resolution Study 

The GEANT detector simulation together with the event reconstruction routine 
allows us to study the effects of secondary scattering in the detector material, energy loss 
and particle mis-identification on the determination of the tracks' kinematics. Figure 17.11 
compares the track transverse momenta and pseudorapidity for the simulated generated 
and the reconstructed primary tracks (electrons and positrons), which were successfully 
associated to the generated tracks. 

From Figure I7~T1 we see that the reconstructed tracks' p r ^ co are distorted compared 
to the true generated values p 9 ± n '. The distortion can be separated into the smearing of the 
reconstructed momentum value for the fixed value of the generated p 9 ^ n and the shift of the 
mean {j ) r ^ co ^ mean compared to the original value. The smearing can be attributed to the 
secondary scattering of the electron/positron tracks in the material of STAR. The shift of 
the mean value is due to the use of the pion mass hypothesis in the fitting of the primary 
tracks and in the corrections for the energy loss applied to the primary track fSection l5.1.ci|) . 
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Reconstructed Transverse Momenta (p™") vs. Generated Transverse Momenta (pj™) Reconstructed Pseudorapidities (n/ e ™) vs. Generated Pseudorapidities (r| 9en ) 




Figure 7.1: Contour plots of transverse momentum and pseudorapidity for the reconstructed 
tracks vs. generated values. Solid lines represent equal generated and reconstructed values. 

The shift and smearing are most significant for p 9 ^ n ~ 95 MeV/c, where the mean of the 
distribution p r f co — p 9 ^ n is equal to 8.5 MeV/c (8.9% of the value of p g ^ n ) and the standard 
deviation of the distribution is 10.3 MeV/c (10.8% of p 9 ± n ). 

The distortion in ry is much less significant. The mean of the distribution r f eco — v p en 
is zero for all values of rf en , and the standard deviation of the distribution is 0.017. For a 
typical track with pseudorapidity of rj = 0.6, this represents only a 2.8% distortion. 

We can correct the data for the shift in the mean value of reconstructed transverse 
momenta. We slice the simulated reconstructed transverse momentum spectrum between 
50 and 150 MeV/c into 20 slices of p r f c ° and find the mean value Sp±^ = (p r f co — P± n ) i 
for each slice. Figure 17.11 shows all 20 values of 8p±i vs. p r f c ° and also a fit to the 4-th 
order polynomial function f(p±). Each value of the data p± is then corrected by the value 
of f(p±)- For a typical value of p± ~ 75 MeV/c, the uncertainty of the correction f(pi_) is 
on the order of 10% of the value of the correction. 

7.2 Efficiency Corrections 

7.2.1 Reconstruction Efficiency 

We would like to find the reconstruction efficiency for e + e~ pairs in the defined 
acceptance region. From the analysis chapter we established the following definition of the 



81 




Pjf co , GeV/c 



Figure 7.2: Value of 5p± j vs. pj 6 ^° (dots) with statistical uncertainties. Solid line - 4-th 
order polynomial fit f{pi_)- 

acceptance region: 

\i] e - \ < 1.15 , p± e - > 65 MeV/c , \r/ e +\ < 1.15 , p ±e + > 65 MeV/c ^ 
140 MeV < M inv < 265 MeV 

For this acceptance region, we will compute the efficiency correction as: 

1 ^generated events in acceptance region 

Reconstruction Efficiency ^events after detector simulation and all analysis cuts 

(7.2) 

This efficiency correction represents the convolution of the efficiencies at several 
stages of analysis, such as tracking, vertex finding and pair identification for the tracks 
in the defined acceptance region. These efficiencies are not independent from each other, 
therefore the reconstruction efficiency cannot be computed as a simple product of separate 
efficiencies. Equation 17.21 also includes the smearing effects. 

Figure 1731 shows the dependence of the reconstruction efficiency on the total trans- 
verse momentum (p t ° t ), invariant mass (Mi nv ), polar angle in the center of mass frame 
(cos(0')) and pair rapidity (1"). Efficiencies if Mi nv and p 1 ^ can be fit with smooth func- 
tions e(Mi nv ) and e(p t ? t ), which we will use as efficiency correction functions. The peaked 
shape of the reconstruction efficiency e{Mi nv ) is a product of raising tracking and vertex 
finding efficiency and dropping pair identification efficiency as a function of M{ nv . The 
reconstruction efficiency in p 1 ^ is flat, except for momenta under 12 MeV/c, where the 
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Reco Efficiency vs. Pair Invariant Mass 
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Figure 7.3: Reconstruction efficiency vs. p^ 1 , M inv , Yp air and cos(6>') . Efficiencies in M inv 
and p 1 ^ are fit with smooth functions e(Mj n „) and e(p F t ) (solid lines), p 1 ^ efficiency is also 
presented for the simulations without p± smearing (dashed histogram). 

efficiency drops, mostly as a result of the p± smearing. As a check procedure, we turned off 
p± smearing in the simulations and the resulting reconstruction efficiency in p 1 ^ (dashed 
histogram in Figure shows a much smaller dip for pairs with p 1 ^ < 12 MeV/c. 



7.2.2 Systematic Errors 

Systematic errors in the reconstruction efficiency determination occur because the 
simulations do not perfectly describe the data. The following subsections discuss the sys- 
tematic errors in the reconstruction efficiency from various sources. 



Tracking 

To check how well the simulations describe the tracking performance of the STAR 
TPC, we compare the ratio of ^hits to the number of possible hits on the track ^possible hits 
in data and Monte Carlo. The number of possible hits was calculated by extrapolating 
a track's geometry over the pad plane. The calculation incorporated both the detector 
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Figure 7.4: Distribution of the hits /^possible hits) and the track fit Xper d.o.f- Dots - data, 
solid - Monte Carlo (scaled). 



geometry and features of the track finding routine. This calculation also correctly handled 
daughters from photon conversions that occurred within the tracking volume of the TPC|52|. 
The distribution of the ratio (hit ratio = N hits /^possible hits) is shown in Figure ITU (left). 

The quality of the track fit in data and Monte Carlo is characterized by the 
Xpcr d.o.f of the track fit. Figure I7.4I (right) shows the distribution of Xper d.o.f in data 
and Monte Carlo. To increase the statistics, events were selected with all analysis cuts in 
the table I6.2I except the particle identification cut. The individual track momenta were 
limited to p < 200 MeV/c. The agreement between Monte Carlo and data is good for both 
distributions. 

In a related analysis ( 20 ) the authors studied the variation in the tracking efficiency- 
corrected p± and T) spectra due to small variations in track quality cuts. The systematic 
error was found to be 6.4% for each reconstructed track, which we accept as a tracking 
efficiency systematic error for the present analysis. 



Vertex Finding 

We use data to produce an estimate of vertex finding efficiency for e + e~ pairs 
and compare the result to the efficiency obtained from Monte Carlo (see Section [6.2.4I and 
Figure |6~d^) . In the data, we plot invariant mass spectra for identified secondary e + e~ pairs 
comprised of any two global tracks (diamonds) and for pairs comprised of two global tracks 
that also were primary tracks (triangles). The first sub-set of pairs is larger than the second 
sub-set, since it includes events with no primary vertex. The ratio of the e + e~ pairs with 
the primary tracks to the total number of pairs gives an estimate of the vertex finding 
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Figure 7.5: Left: invariant mass spectrum from secondary pairs (diamonds - pairs comprised 
of two primary tracks, triangles - pairs comprised of any tracks). Right: comparison of 
vertex finding efficiency from the simulations (crosses) and data (dots). 



efficiency in the data. Due to the low statistics, coarse binning had to be used. 

From Figure I7.5I the general shapes of vertex finding efficiency distributions vs. 
pair invariant mass agree for data and Monte Carlo. The most sizeable difference between 
Monte Carlo efficiency and efficiency from data appears to be between 0.14 GeV < M{ nv < 
0.2 GeV. Analysis of the pf*, cos{9') and Y distributions in the identified secondary e + e~ 
pairs shows significant deviations from Monte Carlo spectra, which are not present in the 
distributions of primary e + e~ pairs. We believe that this is due to the fact that without 
a requirement of a primary vertex, the selected secondary e + e~ pair sample contains a 
large fraction of events which consist of electron or positron tracks that suffered extreme 
energy loss (e.g. hard Bremsstrahlung) in the detector material, or the spiraling tracks 
(Section 15 .2|) . The reconstructed momentum of such tracks differs from the true momentum 
by more than 50%, and therefore these events are unusable for the e + e~ pair spectrum 
determination, and are a background to the 'good' e + e~ events. These events contaminate 
the secondary e + e~ pair sample, causing an under-estimation of the vertex finding efficiency 
from data. The primary vertex requirement is essential in rejecting these backgrounds, and 
must be enforced in the event selection. 

Formally, the ratio of the integrated vertex finding efficiency from Monte Carlo to 
the vertex finding efficiency from the data is 1.17±0.12. This suggests a 17% systematic error 
in vertex finding efficiency from Monte Carlo. To correct for the unaccounted backgrounds 
in the secondary e + e~ pairs sample, we take 50% of 17% (8.5%) as an estimate of the 
systematic error in vertex finding efficiency. 
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Pair Identification 

To estimate systematic uncertainty in pair identification efficiency via specific en- 
ergy loss we compute the systematic error on the value of shifted z^ orrected stmu from the 
correction in Equation (|6.5|) . According the Gaussian error propagation law: 

/ X ^corrected simu \ 



(-.data \ 2 o / I \\2 of f ()rr data \ 2 / rjdata \ 2 \ 

( 5 <Z— >) 2 +(5(^))+(Z— -<Z— >) 2 ^^j +[-^Tu-^z mu ) J 

(7.3) 

where errors 5 (Z^ mu ), 5 (Z data ), 5a s ^ mu and 5o d ? ,ta are the fit errors on the determination 
of the mean and width of the Gaussian fits to the data and simulations peaks in Figure 

Given a single track systematic error fiz% orrected stmu i we systematically shift each 
track's z^ orrected simu in the simulation at first by +SZ W and then by —5Z n and compare 
the resulting particle identification efficiencies. The difference between the two curves was 
found to be at most at the level of 2%; therefore the systematic error in pair identification 
is negligible compared to the statistical error. 

Transverse Momenta Mis-Reconstruction 

To study the systematic error due to shift in reconstructed {p r _^ co ) compared to 
the true value of transverse momenta, we re-compute the reconstruction efficiency without 
applying the shift correction to the simulations. The raw Mi nv distribution (which is com- 
puted with p± correction) is then efficiency-corrected with two different efficiency functions 
- one with account of p± shift, and the other without. Figure ITUl compares the two resulting 
efficiency-corrected Mi nv spectra. 

Figure 17.61 shows that the effect of pj_ distortion is most significant in the last 

and first bins. From Section 17.11 the error of the p± correction is ~ 10% of the value of 

the correction. Therefore, the systematic error on the number of events in the first bin is 

10% • ((91 - 43)/ 43) = 11% and in the last bin it is 10% • ((20 - 5)/20) = 7.5%. However, 

the errors in the first and last bin are strongly anti-correlated. Therefore, for the measure 

1 We scale up each error by the Xper d.o.f °f the fit to account for the general disagreement of the distri- 
bution shape in the data and the fit function. For full explanation of the method see standard reference [fM]. 
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Efficiency Corrected Pair Invariant Mass 
t - p^-corrected efficiency used 
— V - p^-uncorrected efficiency used 
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Figure 7.6: Efficiency-corrected M{ nv data spectrum. Diamonds - p± shift correction was 
not used in efficiency computation. Triangles - efficiency includes p± shift correction. 



of the total cross-section systematic error due to p± distortion we use the total number of 
entries in the two histograms: 



r pi distortion 
™syst 



10% ((177- 146)/146) = 2% 



(7.4) 



Asymmetry 



Asymmetry in the number of reconstructed tracks in the East side of the TPC vs. 
the West side of the TPC has been observed in the high-statistics analyses of the hadronic 
events in STAR|2*U|. Since the simulations assume a perfect symmetry between East /West 
sides of the TPC, a systematic error could be introduced in the efficiency determination. 
To test if there might be any asymmetry in observed e + e~ events we compute the ratio: 

N++ - iV__ N++ + iV__ 14 - 22 



A 



N. 



N 4 



0.15 



(7.5) 



. + iV__ + iV + _ 14 + 22 + 52 

where N ++ , N and are the numbers of events in the identified pool of 52 e + e~ 

events with two tracks with positive longitudinal momenta (p z ), two tracks with negative 
p z and two tracks with p z of opposite charges. 

The ratio A is within 2a ~ 2/v52 from zero, therefore systematic error due to 
asymmetry is not appreciable in comparison with the statistical uncertainty. 



Sensitivity to Generated Spectra 

The effect that p± smearing and particle identification have on the the reconstruc- 
tion efficiency depends in part on the spectra distributions generated by event generator 
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EPA Event Generator Pair Total Transverse Momentum Spectrum 
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Fi gure 7.7: Left: EPA generated p^_ distribution is fit with a smooth function F(p£ ) and 



each event in the sample is weighted by 1/F(p? t ). Right: reconstruction efficiency vs. p 1 ^ 
for flat generated spectrum (solid line) and EPA spectrum (crosses). 



(e.g. (p*|°*)9 eri and Mf™). If the physics model of the event generator spectra is wrong, 
there may be a systematic error introduced in the reconstruction efficiency computation. 

To test the sensitivity of reconstruction efficiency to variations in generated spec- 
tra, we can generate an event sample with, for instance, flat p 1 ^ distribution, pass this event 
sample through the GEANT simulation, event reconstruction and analysis, and determine 
reconstruction efficiency in that way. Alternatively, we can use the same event generator 
event sample that we have been using so far (with pair transverse momentum distributed 
according to the probability distribution /(p^)), but compute the efficiency by the formula: 



Reconstruction Efficiency From Flat p° Spectrum 



# weighted events after detector simulation and all analysis cuts 
^weighted generated events in acceptance 

where the weight for each event is inversely proportional to f{p t ° t ). Incidentally, the effect 

of weighting ensures that the denominator in Equation (|7.6j) is a flat distribution. 

Figure 17771 shows how we get weights for each events (left). Using weighting, we 

compute the reconstruction efficiency from a flat p 1 ^ distribution, and compare it to the 

previously obtained efficiency from EPA Monte Carlo (Figure l7.7| right). The difference 

between the two distributions appears only for low transverse momenta (p 1 " 1 " ~ 10 MeV/c). 

The EPA transverse momentum distribution has a peak at p 1 ^ ~ 5-^10 MeV/c and therefore 

smearing the transverse momentum value upwards in the reconstructed events affects EPA 

distribution more strongly than a flat distribution. This explains why the reconstruction 



efficiency is smaller for the EPA distribution than for the flat distribution. For most values 
of the transverse momenta {p 1 ^ > 20 MeV/c) the reconstruction efficiency is not affected 
by the shape of the generated distribution. 

Results of Varying all Analysis Cuts 

We compare the sensitivity of data and simulations to small variations in cuts 
in table 16.21 by systematically varying each analysis cut and observing the change in the 
number of events which pass the new cuts relative to the number of evens that pass the 
original cuts. All cuts were varied by ±2.5%, ±5.0% and ±10%. 2 The relative event 
number variations were found to be very close (differ on the average by ~ 4.5%) in data 
and simulations, therefore the simulations describe the data response to the analysis cuts 
with excellent accuracy. 

7.3 Backgrounds 

Combinatorial Backgrounds 

Random combinations of tracks in low-multiplicity events (e.g. grazing nuclear 
collisions) can occasionally yield an event which satisfies all analysis cuts. We use the like- 
sign pairs (events with primary pair total charge Qtot 7^ 0) as a model of combinatorial 
backgrounds. Figure shows a total pair transverse momentum spectrum of such pairs. 
The majority of like-sign pairs are above p 1 ^ > 0.15 GeV/c, and there are very few like-sign 
events below p 1 ^ < 0.15 GeV/c, where the coherent e + e~ events are concentrated. After 
applying cuts 0.14 GeV< Mi nv <0.265 and p^ < 0.1 GeV/c to the data, there is one 
background event left. Therefore we neglect the effect of combinatorial backgrounds in the 
data. 

Coherent Pion Background 

A background to ultra-peripheral e + e~ events arises from ultra-peripheral tttt 
pairs, which come from coherent photo-nuclear p° production or direct tttt production (Sec- 
tion 12.4.2(1 and are mis-identified as e + e~ pairs. To estimate the number of 7r + 7r~ pairs 
which might have passed the analysis cuts, we used a Monte Carlo for this process. 12,000 
2 The cut on the total primary pair charge was not varied. 
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Figure 7.8: Total pair transverse momentum and invariant mass of like-sign pairs. Left: 
events selected with 0.14 GeV< M inv <0.265 GeV. Right: events selected with < 0.1 
GeV/c. 



events were generated with 0.12 GeV < Mf n ^ < 0.3 GeV and \Y e e | <1.2. From this 
sample 5 events were found to pass all analysis cuts in the table IB~2l therefore reconstruction 
efficiency for coherent tttt pairs is = 0.0042%. Total expected tttt pair contribution in 
the selected pool of 52 events is then cauAu->Au*Au*+p°+ direct tttt " LjotEff*" 9 • E T7r = 0.026 
events. The value of the cross-section &AuAu^Au*Au*+p° was taken from the theory (5.3 
mb)[Jj and the value of L^Eff*™ 9 was set equal to the total luminosity times trigger effi- 
ciency in our data sample. 

The tttt contribution is extremely low because of a very small fraction of the pairs 
with invariant masses (in e + e~ hypothesis) in the range of 140 to 265 MeV and a strong 
suppression of these pairs by particle identification requirement. 



Other sources of backgrounds 

Coherent ultra-peripheral p° mesons can decay into e + e~ pairs. Due to a com- 
bination of the small p° production cross-section and a small p° — ► e + e~ branching rate 
(~ 5 • 10~ 5 ) this reaction does not produce any appreciable backgrounds. 

Coherent p + p~ are not expected to contribute any backgrounds due to the small 
cross-section of this reaction {a^/a ee ~ (m e /m M ) 4 ). 

Individual protons inside the heavy ions can emit photons independently of the 
rest of the protons in the ion. These photons and their interactions are called 'incoherent', 
to distinguish them from the photons emitted by the ions a a whole. The cross-section for 
incoherent two-photon processes is suppressed by a factor of 1/Z 2 ~ 1.6 • 10 -4 , therefore 
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their contribution is negligible. 

Cosmic muons can be reconstructed as a pair of tracks in STAR TPC with the 
total charge zero. These spurious pairs are strongly suppressed by the requirement of 
the coincident neutron signals in the ZDCs. Additional suppression of this background is 
provided by the vertex cut R ver t < 9 cm and \z ver t\ < 100 cm. The remaining cosmics rate 
is negligible. The vertex requirement also rejects the pairs consisting of the tracks which 
come from the beam-gas interactions. 

Photons produced in hadronic AuAu collisions can pair-convert into an e + e~ pair 
in the material of SVT (located at ~ 6 cm -j- 15 cm and in the STAR TPC gas vol- 
ume). However, the opening angle of the e + e~ from photon conversions is very small 
(~ 2m e /E p hoton ~ 0-01 radians for a photon of energy E p h oton = 100 MeV) and the invari- 
ant mass of such pairs is extremely small, well outside of our acceptance range 65 . 

7.4 Efficiency Corrected Spectra 

We correct the raw event spectra by the reconstruction efficiency Figure I7U1 shows 
the resulting Mi nv , Y, cos(9') and p 1 ^ 1 distributions. 

7.5 Extrapolation to 47r 

The acceptance cuts p*™ cfe > 65 MeV/c and \rf rack \ < 1.15 are specific to the 
STAR detector. We would like to extrapolate the observed cross-section to the full 4tt 
range. The kinematic range of the e + e~ pair (i.e. M{ nv and Y) will need to be limited, 
since we do not want to extrapolate to the regions where we have no detecting power. 

First of all, the electron/positron tracks are required have transverse momenta of 
more than p™ n = 65 MeV/c and the tracks are nearly back-to-back in the transverse plane. 
This means that the invariant mass of the pair should be greater than 2p™" = 130 GeV. 
The analysis cut M{ nv > 140 MeV already ensures that the pair invariant mass is above this 
cut-off. The invariant mass of the pair is also limited from above: M{ nv < 265 MeV/c due 
to dropping pair identification efficiency at high masses. We have to carry this cut on to 
the pair kinematics cuts. The second cut we would like to translate into the pair kinematics 
is the cut \rf rack \ < 1.15. As we have shown in Chapter |3] this limits the total pair rapidity 
to \Y\ < 1.15. This is a final pair kinematics cut. 
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Figure 7.9: Resulting efficiency corrected spectra. 



Below is a summary of the previously used cuts ('detector acceptance cuts') and 
the limits of the pair kinematic range we would like to extrapolate to ('47r kinematic range'). 



STAR Detector Acceptance 


Air Kinematic Range 


7/ e - <1.15 , p ±e - >65 MeV/c , \r] e + 1 < 1.15 , p ±e + >65 MeV/c 
140 MeV< M inv <265 MeV 


|Y|<1.15 

140 MeV< M inv <265 MeV 



Table 7.1: Comparison of the detector acceptance cuts and pair kinematics cuts. 



Monte Carlo distributions show that all events passing the detector acceptance 
group of cuts are within the 4-7T kinematic range. The reverse is not true, however. The 
47r extrapolation factor is the inverse ratio of the number of events passing the detector 
acceptance cuts to the number of events within the 4tt kinematic range. 

Figure [?. 1UI shows the ratio of the number of events in the detector acceptance to 
the number of events in the 4n kinematic range and the 4tt extrapolation factor as functions 
of pair invariant mass. The extrapolation factor is rather large. This is due to the fact that 
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Figure 7.10: 4-7T extrapolation factor vs. pair invariant mass. 



the electron and positron tracks are very forward peaked (see Figure As a result, for 
the pairs in the central rapidity range (within the 4tt kinematic range of |V| < 1.15) a large 
fraction of electron or positron tracks has high values of pseudorapidity (rj ~ 4) outside of 
the detector acceptance range. 

As a function of pair invariant mass, the extrapolation factor is larger for the lower 
invariant masses. By restricting the individual track transverse momentum to p± > 65 MeV 
we suppress pairs with invariant masses of the order of Mi nv ~ 130 -j- 180 MeV, but this 
doesn't affect the pairs with higher invariant masses Mi nv ~ 180 -j- 265 MeV, which are 
typically comprised of higher transverse momentum tracks (p± ~ 100 MeV/c). 



93 



Chapter 8 

Results 

We present the differential cross-section distributions of the e + e~ pairs, which 
we obtain normalizing the e + e~ pair spectra by the total luminosity, according to the 
Equation (|6,lj) . We compare the experimental results to the event generator results, and to 
independent theoretical predictions. 

8.1 Cross-section in Detector Acceptance 

Figure lSTTl shows the measured differential cross-sections da e + e - /dMi nv , da e + e - /dp t 
da e + e - l&Y and da e + e - / dcos{6') and compares them to the cross-sections from the event 
generator ('EPA Monte Carlo') and the QED calcuation. 

The shape of the measured distributions da/dMi nv , da/dY and da/dcos(9') agrees 
well with the EPA Monte Carlo distribution, with all data points lying no more than 2a 
away from the Monte Carlo. The distribution of da e + e - / 'dp 1 ^ also agrees well with the EPA 
Monte Carlo, except in the lowest bin. In that bin, for e + e~ pairs with p 1 ^ <15 MeV/c, 
the measured cross-section is 4.5a away from the Monte Carlo prediction. However, due to 
very low statistics (number of raw events in the first bin is 6), the Gaussian '2<r deviation' 
criterion does not apply, and we cannot conclusively reject the possibility that EPA Monte 
Carlo describes the data accurately. 

In the STAR detector acceptance 

\ Ve - \ < 1.15 , P±e - > 65 MeV/c , |r/ e+ | < 1.15 , p ±e+ > 65 MeV/c 
140 MeV < M inv < 265 MeV 
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STAR eV Differential Cross-section vs. Pair Invariant Mass 
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Figure 8.1: Comparison of measured electron-positron production cross-sections (dots) with 
event generator output (solid histogram), and lowest-order QED calculation (dashed). Error 
bars incorporate statistical and systematic uncertainties. 
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we find the cross-section of ultra-peripheral e + e production with mutual nuclear excitation 
at RHIC to be: 

<TAuAu-+Au*Au*+e+e- = L65 ± 0.23(stat.) ± 0.30(syst.) mb (8.2) 

which is 1.2a less than a EPA Monte Carlo prediction of 2.08 mb. The systematic un- 
certainty quoted consists of a 10% contribution from luminosity normalization, 6.4% con- 
tribution per track from tracking (total 13% tracking systematic uncertainty) and 8.5% 
systematic uncertainty from vertex finding. These are added in quadratures for a total 
systematic error of 18.5%. 

Additionally, we have recently been provided with a theoretical prediction of the 
AuAu — > Au*Au* + e + e~ cross-section at RHIC in the STAR detector acceptance by K. 
Hencken|48j. This computation is based on the full QED lowest-order term for the process 
77 — > e + e~, taking into account photon virtualities and treating the external Coulomb 
fields of the Au ions in the external field approximation ,T3]. The computed differential 
cross-sections are presented in Figure l8~Tl as dashed histograms ('QED calculation'). 

All four differential cross-section distributions in the data show excellent agreement 
with the QED calculation. The e + e~ production cross-section prediction from the QED 
calculation is 1.88 mb, which is only 0.6a higher than the experimental measurement. The 
agreement of the data vs. QED calculation for da e + e - /dp^ distribution is much better than 
the agreement of data vs. EPA Monte Carlo. We believe that this is due to the fact that 
the QED calculation takes photon virtuality and transverse momenta into the account in 
computing the e + e~ production cross-section. The Equivalent Photon Approximation, on 
the other hand, assumes that the photons have zero virtuality and the transverse momenta 
of the photons are ignored in the cross-section computation. 

8.2 47r-extrapolated Cross-section 

Using the extrapolation factor in Figure 17.101 we extrapolated the measured dif- 
ferential cross-section da/dMi nv in the STAR detector acceptance to the differential cross- 
section in the full 4-7T solid angle for the limited kinematic range: 

140 MeV < M inv < 265 MeV , \Y\ < 1.15 (8.3) 
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STAR do e+e . / dM inv : 47i-extrapolated, |Y|<1.15 




Pair Invariant Mass M inv , GeV 

Figure 8.2: 4tt extrapolated cross-section for reaction AuAu — ► Au*Au* + e + e~at STAR. 
Solid histogram - EPA Monte Carlo. 

The cross-section was found to be: 



VAuAu^Au*Au*+e+e- Utt = 40.8 ± 5.7(stat.) ± 7.5(syst.) mb 



U) 



and the EPA Monte Carlo prediction of the e + e~ production cross-section in the 4-7T solid 
angle is 46.4 mb. The extrapolated cross-section and the EPA Monte Carlo prediction are 
shown in Figure |HI2l 
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Chapter 9 

Conclusions 

We observe coherent ultra-peripheral e + e~ pairs with mutual nuclear excitation in 
STAR. The cross-section of e + e~ production agrees with the EPA Monte Carlo simulation 
result and an independent theoretical prediction based on the full QED calculation of the 
e + e~ production. With the available statistics the data is consistent with the assump- 
tion that e + e~ production is independent from the mutual nuclear excitations at RHIC. 
The data is also consistent with the lowest-order two-photon QED approximation of the 
electromagnetic e + e~ production. 

The differential cross-section distributions da/dMi nv , da/dp 1 ^, da/dY and da / dcos(6') 
agree with both the EPA Monte Carlo and QED calculation, except for the low p 1 ^ events. 
For events with p 1 ^ < 15 MeV/c the data favors the QED calculation over the EPA Monte 
Carlo. 

Due to low statistics in the identified e + e~ dataset, we are not able to conclusively 
rule out the validity of EPA approach for e + e~ pair production. A tenfold increase in 
statistics would allow us to answer the question of the EPA applicability at RHIC for e + e~ 
production. Increasing statistics by a factor of 100 would most likely allow us to observe 
the higher-order QED contributions to the e + e~ production cross-section. 

Many improvements can be made in the experimental setup of STAR to enhance 
the collected e + e~ statistics. Lowering the magnetic field to 0.125 T could extend the 
TPC tracking to tracks with low p±, thus greatly increasing the statistics. Forward Time 
Projection Chambers can also provide tracking of low p± or high \rj\ tracks. Using the Multi- 
wire Proportional Chambers for triggering on the high |?7| tracks would provide triggering 
on the e + e~ pairs produced without the mutual nuclear excitations. The cross-section 
of the reaction AuAu — > AuAu + e + e" is much larger than the cross-section of the e + e~ 
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production with mutual nuclear excitation, therefore the statistics for such events will be 
much higher. With the growing interest in the ultra-peripheral physics results within the 
RHIC community, the perspectives of e + e~ production studies in STAR are promising. 
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Appendix A 

Formulae Reference 



Determination of the track momenta from the helix parameters 



Symbols Name (units) Definition 
(xo, yo, zoy coordinate of first point (cm) 

xy-plane direction at first point 

dip angle (°) tan -1 (p z /p±) 
( x cUc) coordinate of helix center (cm) 
Rt helix radius (cm) 

rigidity units R(m) and B(T) 
q* charge (relative to the charge of one e + ) 

p* ± transverse momentum, p xy , (GeV/c) 0.3RqB, units i?(m) and B(T) 

p z momentum along the beam axis (GeV/c) pj_ • tan (A) 

p total momentum (GeV/c) \J p\ + Pi 

E Energy assuming particle of mass(m) (GeV) \J p 2 + m 2 

V rapidity I ln (Sf) 

77 pseudorapidity — In [tan ( 18 °2~ A )] 



Table A.l: Formulae for determination of the track momenta from the helix parameters. 
The symbols f and % denote variables that are obtained directly from the helix parameters 
of tracks and the photon finding algorithm, respectively. The * symbol is used to denote 
that in the tracking algorithm, the magnitude of a particle's charge is assumed to be the 
charge of a positron. Later, a more accurate assumption of a particle's charge can be made 
through the combination of dE/dx and rigidity. 
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Lorenz Transform 



We work in the (1,-1,-1,-1) metric system. If covariant 4-vector p^ is equal 
(E,p x ,p y ,p z ) then contravariant 4-vector p^ is equal (E, —p x , ~Py, ~Pz)- The Lorentz trans- 
form tensor as a function of (3 = ((3 X , (3 y , (3 Z ) is given by: 



hp: 



( 7 -7/3, 

-jf3 x 1 + TJ 



IPy jp 

(7-l)fefe 



P 1 , 
1 + ^ 

( 7 -l)^fe 



1 + (7-1)^ 
1 + ^ 



\ "7/3; 

where (3 = ^ (31 + (3% + 0* and 7 
The Lorentz transform of vector p 11 is: 

Electron-Positron Pair Kinematics 

Let 

Pi = (\Pi\ ,Pi) = (pi,l)i±,Piz) = (pi,Pix,Piy,Piz) and 

P2 = (\P~2\ ,P~2) = (P2,l?2±,P2z) = (p2,P2x,P2y,P2z) 

be electron and positron 4-momenta (correspondingly) in the lab frame 1 . Lorentz transform 
to the center of mass frame is specified by the boost factor 

yCM = _{~Pl+~P2) 



Pi +P2 

In the center of mass frame electron/positron momenta are given by 

Wf M = At (? CM ) - Pi = (P? M ,P? X M ,P% M ,P? Z M ) and 

Wf M = A"„ (l3 CM ) -rf = (p? M ,pg M , P g M ,pg M ) 
Pair kinematical variables can be then expressed as: 



pY 



P 1_L + P 2±\ 



1 Electron/positron mass is ignored compared to the momentum. 
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Y'pair _ _ } n 



(pi +P2) ~ {Plz +P2z 

(Pi +P2) + (Plz +P22 



Afcnw = ypi -p 2ai 

1 / ( "P IX + ~P2±) 



p± 

, -1 / ( ~Pl_L + "P2_0 • ~P 2_L 
W = COS j-t pz; j 

V P±-\P2±\ 

cos" 1 (pg*/pf M ) , ifpg" >0 
180° - cos" 1 (pf 2 M /pf A/ ) , if pg M < 

Euclidean Rotation 

A Euclidean Rotation can be specified by an axis of rotation ~u = (u x , v y ,v^) and 
angle of rotation (j> around the axis v . Then a transformation matrix R is 

^ z/ 2 + (1 — v^) cos 4> v x v y (l — cos 4>) — v z sin <f> v x v z (\ — cos </>) + i/y sin \ 
R(u,cj))= v x v y {l — cos 0) + z/g sin z/ 2 + (1 — i/ 2 ) cos v y v z (l — cos(f>) — i/jSint/) 

\ v x v z {\ — cos<fi) — i/y sin0 v y v z (l — cos 0) + 1^ sin <j> v 2 z + (1 — v 2 z ) cos y 

A rotation of vector "a* = (a x ,a y ,a z ) can be specified as 

^rotated = R .^ 

Effects of Photon Polarization 

Cross-section formulae (j2.5p and (|2.6[) assume that the photons in the center of 
mass frame are circularly polarized. If this is not the case, we need to distinguish two cross- 
sections: scalar cross-sections a s (two photon polarizations are parallel) and pseudoscalar 
cross-section a ps (perpendicular polarizations). For two photons in the lab frame with 
momenta ki = (0, 0, u)\), k 2 = (0, 0,0*2) and polarizations T\ = (ei_L, 0) , = (e%\_, 0) the 
total cross-section is then|15j: 

a(u 1 ,u}2,e 1 ,e 1 ) = (ex-e 1 ) (T s (wi,u 2 ) + (ei x ei) a ps (o>i, u 2 ) (A.l) 



where a s and a ps are given bv|15j: 



2 1 / 4m 2 12m 4 \ / yfs j s 7\ 2 / 1 6m 2 \ / 4m 2 
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a ps _ (s) = 8ira 2 - ( 1- 

77— >e+e w 1 



4m 2 4m 4 



In I — + / S 

2m V 4m 2 



1 +4vra 2 



1 2m 2 



Am 2 



In the range of center of mass energies where STAR can observe e + e pairs, the 
two cross-sections are approximately equal. In this case, the expression (jA.lj) evaluates to 

a (ui,U2, et, ei) = (e? • e^) 2 + (ei x e^) 2 • <r (s) = a (s) 



since polarization vectors have unit length (\~e 
two cases of polarization separately. 



1). Therefore we do not need to consider 



Two-Photon Luminosity 

We start with an expression defining the two-photon luminosity 

dL~ 



& AuAu^AuAu+e+ e~ ~ 
on the other hand, cr AuAu _ fAuAu+e+ 



a 



77— >e~^e~ 



(M inv ,Y)dM inv dY 



(A.2) 



dM inv dY 

can be computed directly from photon densities 1)2.8)) 



n-- I I I I cr 77 (M» nt , (u;i,u^) , Y" (wi, w 2 )) ~=z~ ( b ( b 2,^2) ^-—*-d 2 b 1 d 2 b 2 




cr 77 [u>i, 0J2) 



Ryj (loi,lo 2 ) 
UJ1UJ2 



d? b 



duj\dw2 



U>1 UJ2 



(A.3) 



Using the change of variable theorem, we can convert integral in Equation HA,3|) 
into the integral over variables (Mi nv ,Y) using the conversion (cji,u 2 ) = (Mi nv exp(y)/2, Mi nv exp(— Y)/2) 



a- 



CT 77 (ui (M inv , Y) , uj 2 (M inv , Y))- 



F 77 (m (M inv , Y) , L02 (Mjrw, y)) 
u 1 (M inv ,Y)-u 2 (M inv ,Y) 



d(u)i,u) 2 ) 



where the Jacobian 

8(uJi,UJ2 



(uJl,UJ2) 



d{M inv ,Y) 



is equal to: 



d(M inv ,Y) 



duJi du>i 

dM lnv dY 

du)2 du>2 

dM mv dY 



cxp(y) M inv cxp(y) 

2 2 

cxp(-y) -Minv exp(-y) 

2 2 



d(M mv ,Y) 



dM im ,dY 



Thus, Equation (|A.3|) can be re-written as: 

F 77 (wi (M in „, Y) , uj 2 (M inv , y)) Mi, 



cr AuAu^AuAu+e+ e 



M 2 /4 

invl 



■a. 



77— >e" r e" 



(M ira; ,y)-(2M im) dy 



Comparing this to the Equation (|A.2|) we can conclude 



dh, 



77 



dM im ,dY 



2F 77 (o^i (M mt) , y) , Lo 2 (Minv, y)) 



